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Alexei Borodin and Grigori Olshanski 

Abstract. We introduce and study a family of Markov processes on partitions. 
The processes preserve the so-called z-measures on partitions previously studied in 
connection with harmonic analysis on the infinite symmetric group. We show that the 
dynamical correlation functions of these processes have determinantal structure and 
we explicitly compute their correlation kernels. We also compute the scaling limits 
of the kernels in two different regimes. The limit kernels describe the asymptotic 
behavior of large rows and columns of the corresponding random Young diagrams, 
and the behavior of the Young diagrams near the diagonal. 

Our results show that recently discovered analogy between random partitions aris- 
ing in representation theory and spectra of random matrices extends to the associated 
time-dependent models. 



Introduction 

In a series of papers (see [BOl], [012], references therein, and also [B05]) we 
have been studying a remarkable family of probability distributions on partitions 
(equivalently, Young diagrams) called z-measures. These objects have a represen- 
tation theoretic origin, they arise in harmonic analysis on the infinite symmetric 
group, see [KOVl], [K0V2]. Surprisingly enough, the z-measures turned out to be 
related to a number of probabilistic models of random matrix theory, stochastic 
growth, random tilings, percolation theory, etc. In this paper, we introduce and 
study a family of Markov processes on partitions which preserve the z-mcasures. 
Our main result is the computation of the dynamical correlation functions for these 
Markov processes. We also compute the scaling limits of the correlation functions 
corresponding to two different limit regimes as the size of partitions tends to in- 
finity. In the first regime we look at the largest rows and columns of the random 
Young diagram^ while in the second one we focus on the boundary of the random 
Young diagram near the diagonal. 

Examples of dynamical models of random matrix type are well known. The 
sources of dynamics may be very different: in the Gaussian random matrix ensem- 
bles one allows the matrix elements to evolve according to the stationary Ornstein- 
Uhlenbeck process (Dyson's Brownian motion [Dy]), in tiling models one reads the 
two-dimensional picture section by section [Jo2], [Jo5], [Jo6], [OkR], in growth 
models the time parameter is present from the very beginning [PS]. 

In our setting the construction of dynamics is different; it is based on represen- 
tation theory. We heavily rely on the fact that the z-measures define characters 
of the infinite symmetric group and thus possess a special coherency property. It 
refiects the consistency of restrictions of a character of the infinite symmetric group 



^This limit regime has a representation theoretic meaning, see [B02], [012]. 
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to various finite subgroups. The resulting Markov processes arc analogous to those 
arising in other models, and in degenerations they even coincide with some of those, 
see [BO 7]. It is rather surprising that the similarity of the z-measures to measures 
of different origin extends to dynamics associated with those models. 

One of the elements of our construction is a special family of birth and death pro- 
cesses associated with Meixner orthogonal polynomials. Such birth and death pro- 
cesses, among many others, were extensively studied by Karlin-McGregor [KMGl], 
[KMG2]. Certain degenerations of our Markov processes admit a natural descrip- 
tion in the language of Karlin-McGregor, sec §7.5 below. 

Let us now describe our results in more detail. 

Let Y denote the set of all Young diagrams. Wc consider a family Mz,z',(, of 
probability measures on Y which depend on two complex parameters z and z' and 
a real parameter ^ e (0,1). The weight of a Young diagram A with respect to 
Mz,z',i is given by 

M,,,,,j(A) = (l-0^^'Cl'l {z)x{z') 



dim A 



Here 



wa= n 

(product over the boxes of A) is the generalized Pochhammer symbol, and dim A is 
the dimension of the irreducible representation of the symmetric group of degree 
|A| associated to A. In order for Mz,z'.^{^) to be nonnegative for all A G Y, wc need 
to impose certain restrictions on z and z', for instance, z' = z. All possibilities for 
{z, z') are given before Proposition 1.2 below. 

Let Y„ denote the set of all Yoimg diagrams with n boxes. Restricting Mz,z',( to 

(n) 

Y„ C Y and renormalizing it, we obtain a probability measure on Y„, which 

does not depend on ^. The measure Mz,z',£, may be viewed as a mixture of the 
finite level measures Af^"-*,. 

The Markov processes that we arc about to construct, are jump processes with 
countable state space Y and continuous time t e M. The jumps are of two types: 
one either adds a box to the random Young diagram, or one removes a box from 
the diagram. 

The event of adding or removing a box is governed by a birth and death process 
A^c,{(i) := |'^(*)| on Z+. This process depends on ^ and the product c = zz' , and 
its jump rates are given by 



Prob{7Vc,|(i + dt) = n+l\ Nc,i{t) = n} 
Prob{iVc,5(t + dt) = n-l\ A^c,{(i) = n} 



^{c + n)dt 
ndt 



This is special case of the birth and death processes considered in [KMG2]. Its 
invariant distribution, the so-called negative binomial distribution, is the weight 
function for the Mcixncr orthogonal polynomials. 

Conditioned on the jump n ^ n+ 1, the choice of the box to be added to 
A is made according to the transition probabilities 

p' (n, X;n + l,iy)= ! — , 

(zz' + n)(n + 1) dim A 
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and conditioned on the jiunp n ^ n — 1, the choice of the box to be removed 
from A is made according to the cotransition probabilities 

I dim^ 

The transition and cotransition probabihties are naturaUy associated with finite 
level measures M^^"l . These probabilities were introduced in [VK] in the context of 
general characters of the infinite symmetric group (sec also [Kc2]). 

The jump rates X y v and A \ correctly define a stationary Markov process 
^z,z',i{t) on Y. The measure Mz^z',i is the invariant measure for this process. 
Moreover, ^/^^ is reversible. In the degenerate case of z or z' being an integer, 
^z,z',i can be interpreted in terms of finitely many independent birth and death 
processes subject to a nonintersection condition, see §7.5 below. 

One can also construct Markov chains which preserve the finite level measures 

in) 

. The key idea is that finite level measures are preserved by transition and 
cotransition probabilities. Thus, adding a random box and removing a random 
box afterwards leaves M^^, invariant. Alternatively, one can first remove a box 
and then add a box. These two procedures yield two different Markov chains. 
They were suggested by Kerov a long time ago (unpublished). The same idea was 
independently exploited by Fulman [Fu]. It should be noted that our methods 
based on dctcrminantal point processes arc not directly applicable to such Markov 
chains. The idea of mixing all finite level measures together-^ is essential for us, 
it allows us to obtain explicit formulas for dynamical correlation functions, as we 
explain below. 

It is well known that Young diagrams can be viewed as infinite subsets (point 
configurations) in a one-dimensional lattice. This parametrization of Young dia- 
grams turns out to be very useful. 

Let U be the lattice of (proper) half integers 

7' — 7 A- 1 — i _5_3_i 1 3 5 \ 

'^'1 \- • ■ 1 2' 2' 2' 2' 2' 2''"'J'" 

For any A e Y we set 

X(A) = {Ai-i + i |i = l,2,...}cZ'. 

For instance, for the empty diagram A = 0,X(0) = {...,— Using the 
correspondence A i— > X(A) we interpret the measure M^^z'^i on Y as a probability 
measure on 2^ . This makes it possible to speak about the dynamical correlation 
functions of Az^z',^ which uniquely determine the process. They are defined by 

Pn{ti,xi;t2, X2;- ■ tn, Xn) = Prob {X_ (A) at time ti contains for 1 < i < n} . 

Here n = 1,2,..., and the nth correlation function p„ is a function of n pairwise 
distinct arguments (ti, xi), . . . (t„, a;„) G R x Z'. 

The notion of the dynamical correlation functions is a hybrid of the finite- 
dimensional distributions of a stochastic process and standard correlation functions 
of probability measures on point configurations. 



^which may be viewed as a passage to the grand canonical ensemble, cf. [Ve] 
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The reason why we are interested in dynamical correlation functions is the same 
as in the "static" (fixed time) case: As we take scaling limits of our processes, 
the notion of weight of a point configuration ceases to make any sense because the 
space of relevant point configurations becomes uncountable. On the other hand, 
the scaling limits of the correlation functions do exist, and they carry complete 
information about the asymptotic behavior of our processes. 

Theorem A (Part 1). The dynamical correlation Junctions of ^z,z',i, have the 
determinantal form (n = 1,2,...) 

Pnitl . . .] Xn, tn) = dct z,z' ,^ iti-, Xi \ tj , a^j)] ^ , 

where the correlation kernel K_z,z',i{s,x;t,y) is a function on (M x Z')^ which can 
be explicitly computed. 

One way of writing the kernel is by a double contour integral 

Kz,z',ds,x;t,y) 

ei(^-t)T{-z' -x+ ^)T{-z ~y+ \){-ir+y+^ 
(r(-z -x + \)n-z' - X + \)V{-z -y+ \)T{-z' -y+\)Y 

{uJl} {<^2} 

UJl^ '^1^2^ ^ d,UJiduJ2 

e«-* (^'1 - ve) (-'2 - v^) - (1 - v^^i) (1 - ve^2) 

with the contours {oJi} and {us^} o/wi and u>2 satisfying the following conditions: 

• {oJi} and {u)^} go around in positive direction and pass between \/f and 

• The contours are chosen so that the denominator in the formula above does 
not vanish. There are two possibilities of doing that; one of them is used for s>t, 
and the other one is used for s <t, see Theorem 7. 1 below for details. 

This integral representation is convenient for computing the scaling limits of the 
correlation functions. However, it does not reveal important structural features of 
the kernel. Let us now present another way of writing the correlation kernel. 

Consider a second order difference operator on Z', depending on parameters 
{z,z' , S.) and acting on functions /( • ) G i?2(Z') as follows 

iDf){x) = ^C(^ + a; + i)(^' + a:+i)/(x + 1) 
+ ^az + x-^){z' + x-^)f{x -l)-{x + ^{z + z' + x)) f{x). 

This is a self-adjoint operator with discrete simple spectrum SpD = {(1 — 0^'}- 
Its eigenfunctions ■i/'o) 

D'4}a = (1 - 0« • V'a, 

are explicitly written through the Gauss hypergeometric function, see (2.1) below. 
We normalize them by the condition HV'all = 1- 
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Theorem A (Part 2). The correlation kernel for the dynamical correlation func- 
tions of the Markov process i^z,z',(, can also he written as 

Kz,z,,i{s,x;t,y)=± ^ e-»l«-*lv±a(a;) V±a(y) 

„— 1 3 5 
"~ 2 > 2 > 2 >••• 

with "+" taken for s >t and "— " taken for s <t. 

The functions {ipa} form an orthonormal basis in ^^(Z'). Tims, for s ^ t the 
kernel K_z,z',^ defines a projection operator whose range is the span of the eigen- 
functions of D corresponding to the positive part of the spectrum of D. In this case 
(sec Comments at the end of §3) the kernel can be written in a simpler, so-called 
integrable form: 

P{x)Q{y) - Q{x)P{y) 

= ' 

where P and Q are expressed through the Gauss hypcrgcomctric function. 

The formula of Theorem A (Part 2) shows that our Markov process is determined 
by the following data: a state space 3£, a Hilbert space H of functions on X, a self- 
adjoint operator D in H , and two complementary spectral projection operators P± 
for D. In our case, X = Z', iJ = ^^(Z'), D is the difference operator given above, 
and P± are projections on the positive and negative parts of the spectrum of D. 

It seems that generating Markov processes with dcterminantal correlation func- 
tions by data {X, H, D, P±) of this type is a rather general phenomenon. Similar 
structures have appeared earlier in the dynamics arising in polynuclear growth mod- 
els [PS], [Jo3], in tiling models [Jo5], [Jo6], and in random matrix theory [NF], [Jo4], 
[TW] . Following the terminology of those papers, we call the kernel of Theorem A 
the extended hypergeometric kernel. 

The reader might notice that in our Theorem A as well as in all the papers 
cited above, the values of the extended (dynamical) kernels are always given by 
somewhat different expressions depending on the relative order of the time variables. 
This dichotomy is unavoidable because of a certain discontinuity of the dynamical 
correlation functions. For example, we must have 



s Rj t p2{s,x;t,y) 



P2{s,x;s,y), Xy^y, 
pi{s,x), x^y. 



If we assume the dcterminantal structure of the dynamical correlation functions 
with a kernel K{s, x; t, y) then 



pi{s,x) = K{s,x;s,x), p2{s,x;t,y) = 



K{s,x;s,x) K{s,x;t,y) 
K{t,y;s,x) K{t,y;t,y) 



If we further assume that the kernel is continuous in s, i subject to the condition 
s > t,^ then the above relations imply 

K{s,x]s,y) = lim K{s,x; s - e,y) = lim K{s,x] s + e,y) + 6xy 

e— >+0 e— >+0 



^We could have used s <t equally well, this is a question of convention. For instance, trans- 
position of the kernel does not affect the correlation functions, and this operation turns s >t into 
s<t. 
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The validity of the last relation for the extended hypcrgcornctric kernel can be 
immediately observed from Part 2 of Theorem A using the fact that {ipa} form an 
orthonormal basis. 

Let us now describe our results on scaling limits of the dynamical correlation 
functions. In our previous works we considered three asymptotic regimes for random 
Young diagrams without dynamics: one for largest rows and columns, one for rows 
and columns of intermediate growth, and one for the behavior of the boundary of 
the Young diagrams near the diagonal, see [B05] and references therein. In all 
three limit regimes the parameter ^ tends to 1, which makes the expected number 
of boxes in the random Young diagram go to infinity. 

In this paper we concentrate on the first and the third limit regime, but with the 
presence of dynamics. Let us start with the behavior of large rows and columns. 

In order to catch the largest rows and columns in the limit ^ y 1, we need to 
scale them by (1 — ^). This leads to scaling of the state space Z' by the same factor. 
That is, Z' is replaced by (1 - ^)Z' which in the limit turns into R* = R \ {0}. 

The paramctrization of Young diagrams by point configurations X (A) is not 
suitable for this limit transition. Or, rather to say, the positive part of X (A) indeed 
reflects the behavior of largest rows, while the behavior of the largest columns is 
captured by the complement of the negative part of ^ (A) in {...,— 
Thus, instead of encoding A by X (A) we use the map 

A^X(A)=(x(A)n{i,i,|,...})u({...,-|,-i,-i}\X(A)). 

We refer the reader to [B02], [012] for representation theoretic interpretation of 
this map and for further details. 

Theorem B. The scaling limits, as ^ — > 1, of the dynamical correlation functions 
of ^z,z',(, corresponding to the map A i-^ X{X), under the reseating of 11 by (1 — 
have determinantal form with the correlation kernel ii"]^, on (M x M*)^. 

This kernel has four blocks according to the choices of signs of u and v. 
The block with u,v > has an integral representation 

K^,,{s,u;t,v) = e^'('-+''\u/v)^ei^'-'^ 

^ ^ cr(i+Circ2~^(i+C2r ^ ^ 



(27ri)2 7 ^ ' ' ' e--*(l+Ci)(l + C2)-CiC2 

+ 00 +00 

with different choices of contours for s>t and s <t, see Theorem 9.4 below. 
The same block has a series representation 

K^,,{s,u;t,v) = ± Yl e-»l^-*l«;±„(«)w;±„(t;), 

"~ 2 ' 2 ' 2 

where "+ " is taken for s>t, "— " is taken for s <t, and 

= lim(l-0^a([(l-0-'«]) 

are eigenf unctions of a second order differential operator on 
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which are explicitly written through the Whittaker functions, see (9.1) below. 

Similar expressions are available for three other blocks of K^^,{s,u;t,v), see 
Theorems 9.2 and 9.4 below. 

We call KYz'{s,u;t,v) the extended Whittaker kernel. 

In the "static" case s = t the kernel admits a simpler "integrable" form, see 
[BOl], [B02], [Bl], [012], and (9.4) below. 

Let us now proceed to the other limit regime which describes the behavior of the 
Young diagrams near the diagonal. This just means that we stay on the lattice Z'. 
For this asymptotic regime it does not really matter whether we use X_ (A) or X{X) 
to encode the Young diagrams. We refer to [BOS] for a detailed discussion of this 
regime. 

In the following statement we will use a more detailed notation Wa{u; z, z') for 
the functions Wa{u) introduced above. 

Theorem C. The limits, as ^ — > 1, of the dynamical correlation functions ofAz^z',^ 
corresponding to the map A ^(A), under the resettling of time by (1 — ^)""'^, have 
determinantal form with the correlation kernel K ?'^™'"'^(g. x\ r. li) on (M x Z')^. 

For a > T, the correlation kernel can be written in two different ways: as a 
double contour integral 

Kf^-'^a^x.T^y) 

T{-z' -x+ \)T{-z -y+ i)e-'^'(^+^')(_i)^+y+i 
(r(-^ - X + ^)T{-z' -x+ ^)T{-z -y+ i)r(-^' - 2/ + i)) ^ 

+00 +00 

and as a single integral 

r+oo 

Kl—{c7, x; T,y)= / e-^^^-^^w^u; -z, -z')wy{u; -z, -z')du. 
Jo 

The values of the kernel for a < t are obtained from the above formulas using the 
symmetry property 

(a,a;;T,y) = (-1)^+^+^^™;(t, -a;; a, -y), a^r. 

For a = T the kernel admits a simpler expression of "integrable" type 

P{x)Q{y) - Q{x)P{y) 
x-y 

where P and Q are are expressed through gamma functions only, see [B05] and 
(10.3) below. That kernel was called the gamma kernel, and for this reason we call 
iff^f ™^(o", .x; r, y) the extended gamma kernel. 

Note that the extended gamma kernel fits into the same abstract scheme as 
the extended hypergeometric kernel: one takes 3i = Z', H = £^{Z'), D is the 
special case of the difference operator given above corresponding to the limit value 
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^ = 1. The spectrum of this operator fills the whole real axis, the eigenf unctions 
are x i— > Wx{u\ —z, —z'), and the spectral projections P± again correspond to the 
positive and negative parts of the spectrum. 

The functions ipa{x) = ^pa{x; z, z' used in the discussion of the extended 
hypergeometric kernel have the following symmetry: 

ipaix; z, z', = il>x{a; -z, -z', 0- 

This means, in particular, that ipa{x) satisfies second order difference equations 
both in a and x (the bispectrality property, see [Gr]). The two limit transitions 
considered above (Theorems B and C) correspond to taking continuous limits in x 
and a, respectively. This explains why we end up with the same functions Wa(u) 
in Theorems B and C. 

Let us make some remarks about our proof of Theorem A. As a matter of fact, we 
prove the theorem in a greater generality. We introduce certain time inhomogeneous 
Markov processes on partitions. Their fixed time distributions are also the measures 
but now ^ = ^(t) varies with time t. The construction of these processes 
is similar to the stationary ones except that the birth and death process on Z_|_ 
becomes time inhomogeneous. In particular, we consider pure birth and pure death 
processes for which the Young diagrams either always gain new boxes or always 
lose their boxes. These "pure" processes are simpler, their transition probabilities 
can be evaluated explicitly. They can also be viewed as building blocks of general 
processes, more exactly, the transition matrix P{s,t) for a general process can be 
represented as a product P{s,t) = P^{s,u)P^ {u,t) of transition matrices of "pure" 
processes for a suitable intermediate time moment u G (s,t). 

This product representation of the transition matrix plays an important role in 
the proof of Theorem A. We first prove the theorem for a degenerate case, when one 
of the parameters z, z' is an integer, and the process is "finite-dimensional" , that is, 
it lives on the Young diagrams with bounded number of rows or columns. Then the 
needed formulas are derived from a version of Eynard-Mehta theorem on spectral 
correlations of coupled random matrices [EM].** The passage from the degenerate 
case to the general one is based on analytic continuation in the parameters z and 
z'. This passage is not trivial since we need to extrapolate from the integer points 
to a complex domain. The needed analytic properties of the dynamical correlation 
functions are derived from the product formula for the transition matrix P{s,t) 
mentioned above. Let us also emphasize that in our approach, the introduction of 
time inhomogeneous processes is necessary for handling the stationary case. 

There is one more subtle issue that we would like to mention here. Generally 
speaking, even for birth and death processes, jump rates do not determine the 
transition matrix uniquely, see e.g., [Fe2, ch. XVII, §10]. Since we want to define 
our processes by their jump rates, we need to ensure the uniqueness. We were 
unable to find suitable results in the literature and, therefore, we were forced to 
invent a special sufficiency condition which was suitable for our purposes, see §4. 

Let us point out that there exists another way of obtaining the dynamical corre- 
lation functions of Theorem A, based on the formalism of infinite wedge Fock space. 
In [Ok2] Okounkov gave an elegant derivation of static (s = t) correlation functions 
(initially computed in [B02]) using a representation of SL{2) by the so-called Kerov 



''Other proofs of this theorem can be found in [NF], [Jo3], [TW], [BR]. 



MARKOV PROCESSES ON PARTITIONS 



9 



operators. We can extend Okounkov's approach to derive the formula of Theorem 
A. This alternative path bears some similarity to the formalism of Schur processes 
of [OkR], [Ok3]. However, the Schur processes seem to be not applicable in our sit- 
uation. Note also that despite the beauty of Okounkov's idea, a rigorous realization 
of this approach would have to overcome certain nontrivial technical difficulties. 

One more important subject that we do not touch upon in this paper, is a 
family of Markov processes on partitions related to Plancherel measures. In the 
limit z, z' ^ oo, ^ ^ 0, zz'^ 9 > 0, the measures Mz,z',i tend to the so-called 
poissonized Plancherel measure on Y with Poisson parameter 6. This connection 
was used in [BOO] to study the asymptotics of the Plancherel measures. Using the 
general scheme presented in this paper, one constructs Markov processes on Y which 
preserve the poissonized Plancherel measures. These processes may be viewed as 
degenerations of the processes considered in this paper. They are equivalent to the 
droplet model of polynuclear growth. Our results on this other family of Markov 
processes and their scaling limits are presented in [B07]. Let us note that the 
analog of Theorem A for those processes can be obtained either by limit transition 
from Theorem A or by using the Schur process of [OkR]. 

The present paper is organized as follows. In Section 1 we introduce the z- 
measurcs, the associated transition and cotransition probabilities, and other no- 
tions related to the Young graph. In Section 2 we study the eigenfunctions Va of 
the second order difference operator D on Z'. In Section 3 we prove the static 
variant of Theorem A using the method of analytic continuation and reduction to 
the degenerate case of integral parameters. In Section 4 we introduce time ho- 
mogeneous and inhomogeneous Markov processes on Y, prove their existence and 
uniqueness, and compute the transition probabilities for "pure" ascending and de- 
scending processes. In Section 5 we evaluate the transition matrices for integral 
values of parameters. In Section 6 we study the analytic nature of the dependence 
of the dynamical correlation functions on the parameters. In Section 7 we prove 
Theorem A first in the degenerate case using Eynard-Mehta theorem and Meixner 
polynomials, and then in the general case using analytic continuation. In Section 8 
we derive the dynamical correlation fimctions of A^.^/.^ corresponding to the map 
A 1-^ AT (A) (as opposed to the map A A. (A) used in Theorem A). In Section 9 we 
prove Theorem B, and in Section 10 we prove Theorem C. 
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1. Z-MEASURES 

As in Macdonald [Ma] we identify partitions and Young diagrams. By Y„ we 
denote the set of partitions of a natural number n, or equivalently, the set of Young 
diagrams with n boxes. By Y we denote the set of all Young diagrams, that is, 
the disjoint union of the finite sets Y„, where n = 0, 1,2, . . . (by convention, Yq 
consists of a single element, the empty diagram 0). Given A e Y, let |A| denote the 
number of boxes of A (so that A € let £{\) be the number of nonzero rows in 

A (the length of the partition), and let A' denote the transposed diagram. 

For two Young diagrams A and /x we write ^ y \ (equivalently, A \ /u) if ;U C A 
and \^\ = |A| — 1, or, in other words, /x is obtained from A by removing one box. 
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The Young graph is tlic graph whose vertices are the elenients of Y and the edges 
join aU pairs A) such that /i A. The Young graph will also be denoted by Y. 
Clearly, ji X implies [j! / X! , so that the transposition operation A A' induces 
an involutive automorphism of the Young graph. 

For any A € Y„, standard Young tableaux of shape A can be viewed as paths 

/ A(^) / ■■■ / A^"^ = A 

in Y. Let dim A be the number of all such paths. A convenient explicit formula for 
dim A is 

dimA= — - [] (Ai-i-Aj +j), A e Y„, 

where N is an arbitrary integer > £(A) (the above expression is stable in N). 
For A e Y„, e Y„_i set 

( dim^ 

p^(n,A;n-l,M)= <^ dimA' ' (1.1) 

L 0, otherwise 

and note that 

The numbers p^{n, A; n — 1, /x) are called the cotransition probabilities of the Young 
graph. 

A family {M^")} of probability measures M^") on Y„, n = 0, 1, 2, ... , is called 
a coherent system on Y if the measures are consistent with the cotransition proba- 
bilities in the following sense: 

^^""'Hm) = E MeY„_i, n = 1,2,.... (1.2) 

AeY„ 

This concept has an important representation theoretic meaning. Namely, there 

is a 1-1 correspondence between coherent systems on Y and normalized positive 
definite class functions on the infinite symmetric group, see [VK], [Ke2], [012]. 

Note that the cotransition probabilities are invariant under the involution A i— > A' 
of the Young graph. Consequently, the push-forward of a coherent system under 

this involution is again a coherent system. 

Example 1.1. The Plancherel measures defined by 

(„) _ (dimA)^ 

form a coherent family of probability measures, see [VK]. 
Let T be the set of all infinite paths in Y of the form 

y a(i) y a(2) /'•••/' a(") /■..., a("' e y„. 
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This is a compact topological space (a closed subset of the product space YI^^q Y„). 

A probability measure on T is called central if for any n = 1,2,... and 
any A e Y„, the mass of each cylinder set consisting of all paths with fixed 
X'^^K .... A^") = A depends on A only (and does not depend on X^^\ . . . , A^"^^)). 

Any coherent system {M^")} generates a central measure M. on T. By definition, 
the mass of the cylinder set mentioned above equals M^"\X)/ dim A. The relation 
(1.2) ensures that A4 is correctly defined. This defines a one-to-one correspondence 
between coherent systems {M*^")} on Y and central measure M. on T, see [VK], 
[Ke2], [012]. 

For any central measure A4, 

p^{n, X;n-l,n) = ProbjA*"-^) = /z | A^") = A}, 

which is a justification of the term "cotransition probability" . 
Assuming M(I-^I)(A) > for all A e Y, set 

p^n, X;n + l,iy)= Prob{A("+^) = u \ A*"' = A}, n = |A|. 

In contrast to p^{n,X;n — 1,/Lt), these numbers depend on M. We call them the 
transition probabilities of the central measure M (or of the corresponding coherent 

system {Af*^"^}). The transition probabilities define A4 and {Af*^"^} uniquely. 
Note an important relation between the transition and cotransition probabilities: 

M(") (A)pT (n. A; n + 1, zv) = pi (n + 1, 1/; n, A)M("+i) {p). (1.3) 

It implies, in particular, that 

r M("+^)(iy) dim A 
p^n,X;n + l,u) = < M(")(A)dimz/ ' ^ ^' (1.4) 
I, 0, otherwise. 

If (1^1) (A) vanishes for some A G Y then the definition has to be slightly mod- 
ified. Namely, let suppAl be the set of those A e Y for which M(I^I)(A) > 0. 
Equivalently, A e supp A4 if the set of paths passing through A has positive mass 
with respect to Ai. Note that A G suppAl implies ^ G suppA^ for all /i y X. 
The set suppTW spans a subgraph of Y (which may be called the support of M), 
and the transition probabilities are correctly defined on this subgraph by the same 
formula (1.3). Again, the initial central measure M. is uniquely determined by its 
support and the transition probabilities. 

Note two useful equations 

M("-i)(/Lt) = ^M(")(A)pH",A;n- l,Ai), (1.5) 

A 

M("+i)(!/) = ^M(")(A)pT(n,A;n+l,i/). (1.6) 

A 

We shall need the generalized Pochhammer symbol {z)\: 

e.(\) 

{z)x^Y[{z-i + 1)k, zeC, AgY, 
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where 

is the conventional Pochhammer symbol. Note that 

wa= n i^+j-') 

(i,j)eA 

(product over the boxes of A), which implies at once the symmetry relation 

For two complex parameters z, z' set 

Mi3(A) = i^i^, n = 0,l,..., A.Y„, (1.7) 

where dim A was defined in the beginning of the section. The expression (1.7) 
makes sense if {zz')n does not vanish, i.e., if zz' ^ {0,-1,-2,...}. Obviously, 
(1.7) is symmetric in z, z'. 
Note that (see Example 1.1) 

J'^^<^'(>^^=^pLcHerelW- (1-8) 

Let us say that two nonzero complex numbers z, z' form an admissible pair of 
parameters if one of the following three conditions holds: 

• The numbers z, z' arc not real and arc conjugate to each other. 

• Both z, z' are real and are contained in the same open interval of the form 
(m, m + 1), where m e Z. 

• One of the numbers z, z' (say, z) is a nonzero integer while z' has the same 
sign and, moreover, \z'\ > \z\ — 1. 

Proposition 1.2. If {z,z') is an admissible pair of parameters then {M^^^,} is a 
coherent family of probability measures. 

Proof. It is readily checked that if (and only if) one of the conditions above holds 
then {z)x{z')\ > for all A, sec [B05, Proposition 1.8]. Moreover, (zz')« > 
for all n. Hence (1.7) is nonnegative. The fact that each M^"l is a probability 
measure and the coherency property can be proved in several ways. See, e.g., [Oil], 
[B03]. □ 

We call the measures M^"J the z^m,easures on the floors Y„ of the Young graph. 
Depending on which of the three conditions of Proposition 1.1 holds we will speak 
about the principal, complementary or degenerate series of z-measures, respectively. 
By virtue of (1.8), the z-measures may be viewed as a deformation of the Plancherel 
measure (for any fixed n). The principal series of z-mcas\ircs first appeared in 
[KOVl], see also [K0V2]. For more information about the z-measures and their 
generalizations, see [B02], [BOS], [B04], [BOS], [B06], [Kel]. 

Note that the involution A A' of the Young graph takes M^"^, to 
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Let Aiz,z' be the central measure corresponding to the coherent family {M^J^^,}n=o,i,. 
In the case of the principal or complementary series the support of M.z,z' is the 
whole Y. For the degenerate series it is a proper subset of Y: if ^; = fc = 1, 2, . . . 
and z' > k — 1 then supp Mz,z' consists of diagrams with at most k rows, and if 
z = —k = —1, —2, . . . and z' < — (fc — 1) then supp Adz, z' consists of diagrams with 
at most k columns. 

The transition probabihties of the z-measures are given by 

ri,,(„.A:„ + i..)^ <- + °Wy'fVn^'">'°''' - 

(^r + n)(n + 1) dim A 

where c{v/X) denotes the content of the box = vj\^ that is, c~j — i. Indeed, 
(1.9) follows immediately from (1.4) and (1.7). Note that if A is in supp A^^,^' while 
V is not (which may happen for the degenerate series) then (1.9) vanishes due to 
vanishing of one of the factors z + c{v/\), z' + c{v/\). 

For the Plancherel measure, the transition probabilities are 

PllancHerel (n, A; n + 1, .) = (^^^yy^ , A / 

see [VK]. 

Consider a special case of the negative binomial distribution on Z_|_ depending 
on two parameters a > and ^ £ (0, 1): 

7r„,5(n) = (l-0''^^, n = 0,l,2,... (1.10) 

The next formula defines a probability measure on Y which is the mixture of all 
z-measures M_^"^, with given fixed parameters z, z' and varying n by means of the 
distribution (1.10) on n's, with parameters a = zz' and ^: 

Mz,z',d^)=M^J^'J\X)nzz',m) = {^-0'''^^'Hz)x{z')x (^)' • (l-H) 

We call (1.11) the mixed z-measure. An interpretation of formula (1.11) is given 
in [B05, Definition 1.4]. 

Likewise, consider a mixture of the Plancherel measures, depending on a param- 
eter 6 >0: 

Mp,_.e.e,.(A) =Mi],^l,_,(A)e-^^ = e-^^l^l (^)' . (1.12) 

We call (1.12) the poissonized Plancherel measure. Note that it can be obtained as 
a limit case of the mixed z-measures: 

lim Mz,z',iW = Mpiancherel,eW- 

The main objects of this paper are the z-measures and related Markov processes. 
One can also develop a parallel theory associated with the Plancherel measure. We 
do not pursue this goal in the present paper. An interested reader can found the 
statements of the main results related to the Plancherel measure in our paper [B07] . 
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2. A BASIS IN THE SPACE ON THE LATTICE AND THE MEIXNER POLYNOMIALS 

In this section we examine a nice orthonormal basis in the £^ space on the 1— 
dimensional lattice. The elements of this basis are eigenfunctions of a second order 
difference operator. They can be obtained from the classical Meixner polynomials 
via analytic continuation with respect to parameters. 

Throughout the section we will assume (unless otherwise stated) that {z, z') is in 
the principal series or in the complementary series but not in the degenerate series. 
In particular, z, z' are not integers. 

Consider the lattice of (proper) half-integers 

— W 4- 1 — I _5_3_i i 3 5 \ 

— iLi -r 2 — \- ■ ■ J 2' 2' 2'2'2'2'"""J'" 

We introduce a family of functions on Z' depending on a parameter a G Z' and also 
on our parameters z,z',^: 



(x-z z' c) r^^ + ' + ^^'^^^ + '' + ^^ Vc-^-+"Vi 

''^'''''''^'-\T{z-a+^)T{z'-a+^)) ^ ^ ^' 
F{-z + a+\,-z' + a+\;x + a+l\ ^) 



T{x + a + 1) 



(2.1) 



where F{A, B; C; w) is the Gauss hypergeometric function. 

Let us explain why this expression makes sense. Since, by convention, parameters 
z, z' do not take integral values, r{x + z+^) and r(x + z' +^) have no singularities 
for X G Z'. Moreover, the admissibilty assumptions on {z,z') (see §1) imply that 

T{x + z + ^)r{x + z' + i) > 0, r{z-a+ i)r(z' - a + i) > 0, 

so that we can take the positive value of the square root in (2.1). Next, since ^ € 
(0, 1), we have ^/ (^ — 1) < 0, and as is well known, the function w F{A, B; C\ w) 
is well defined on the negative semi-axis w < 0. Finally, although F{A, B; C, w) is 
not defined at C = 0, -1, -2, . . . , the ratio F{A, B; C, w)/T{C) is well defined for 
all C e C. 

Note also that the functions tpaix'jZjz',^) are real-valued. Their origin will be 
explained below. 

Further, we introduce a second order difference operator D{z,z',^) on the lattice 
Z', depending on parameters z, z', ^ and acting on functions f{x) (where x ranges 
over Z') as follows 

D(z,z',0./(.^) = ^az + x+^){z' + x + ^)f{x + l) 



+ yj^iz + x-^)iz' + x-^)fix -l)-ix + + z' + x)) f{x). 
Note that D{z,z',^) is a symmetric operator in £'^{Z'). 

Proposition 2.1. The functions tpa{x; z, z', ^), where a ranges over Z', are eigen- 
functions of the operator D{z,z',^), 

D{z, z', z, z', = a(l - OV'a(a;; z, z', 0- (2.2) 
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Proof. This equation can be verified using the relation 

w{C -A){C- B)F{A, B; C + 1; w) - (1 - w)C{C - 1)F{A, B;C-1; w) 
+C[C -1-{2C-A-B- l)w]F{A, B; C;w)=0 
for the Gauss hypergeometric function, see, e.g., [Er, 2.8 (45)]. □ 

The next lemma provides us a convenient integral representation for functions 

Lemma 2.2. For any A, B e C, M e Z, and ^ e (0, 1) we have 

F(A,B;M + 1;^) _ r(-A + l)r'^/'(l - 0^ 
r(M + l) ~ T{-A + M+l) 



1 - — uj- 



Here ^ e (0, 1) and {u} is an arbitrary simple contour which goes around the points 
and \/f in the positive direction leaving outside. 

Comments. 1. The branch of the function (1 — ^/Ja;)"^"^ is specified by the con- 
vention that the argument of 1 — \fS,uj equals for real negative values of and 

the same convention is used for the function ^ 

2. Like the Euler integral formula, formula (2.3) does not make evident the 
symmetry A B. 

3. The right-hand side of formula (2.3) makes sense for A — 1,2,..., when 
V[—A + 1) has a singularity. Then the whole expression can be understood, e.g., 
as the limit value as A approaches one of the points 1,2, .... 

Proof. Since both sides of (2.3) are real analytic functions of ^ we may assume that 
^ is small enough. Then we may apply the binomial formula which gives 



-B 



OO 'DO 



^ + l)fc(-B); ^(fe_,_;_M)/2 ^k-l-M 



k\l\ 

k=0 1=0 



After integration only the terms with k = 1 + M survive. It follows that the right- 
hand side of (2.3) is equal to 

(1-0^ ^ T{-A + M + l + l){B)i 



r(-A + M + l) V ^ T(l + M + l)l\ 

^ ' ;>max(0,-M) ^ ' 



We may replace the inequality I > max(0, —M) simply by ? > because for negative 
integral values of M (when we have to start summation from I = —M), the terms 
with I = 0, . . . ,—M — 1 automatically vanish due to the factor r(/ + M + 1) in the 
denominator. Consequently, our expression is equal to 

(1 - O'^Fj-A + 1 + M,B;M +1;^) _ F{A, B; M + 1; ^) 
T{M + 1) ~ r(M + 1) ' 

where we used [Er, 2.9 (4)]. □ 
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Proposition 2.3. We have the following integral representations 



r(x' + ,-. + i)r(.r + ,:' + i)y r{z'-a + ^] 



r{z-a+^)riz'-a+^)J T{z' + x+^) 



{a,} 

and 

1pa{x; Z, z', Oi^aiy; Z, z', = if^^^, (x, y) 



(l - V^a;2)"^"""^ f 1 - .r^-»c.-- (2.5) 

V / V W2 / wi a;2 



where 



Jt{x + Z + l)T{x + z' + i)r(2/ + z + i)r(y + z' + ^) 

^-^'^"'^^ = r(x + .' + i)r(, + . + i) ^'-'^ 

Here each contour is an arbitrary simple loop, oriented in positive direction, sur- 
rounding the points and ^/^, and leaving 1 / outside. We also use the convention 
about the choice of argument as in Comment 1 to Lemma 2.2. 

Proof. Indeed, (2.4) iimricdiatcly follows from (2.1) and (2.3). To prove (2.5) wc 
multiply out the integral representation (2.4) for the first function and the same 
representation for the second function, but with z and z' interchanged. The trans- 
position z ^ z' in (2.4) is justified by the fact the initial formula (2.1) is symmetric 
with respect to z z'. As a result of this trick the gamma prefactors involving a 
are completely cancelled out, and we obtain (2.5) □ 

Proposition 2.4. The functions ipa = tpaix; z, z', where a ranges over Z' , form 
an orthonormal basis in the Hilbert space ^^(Z'). 

Proof. From (2.4) it is not difficult to see that the function il^aix; z, z' has ex- 
ponential decay as a; ^ ±oo. In particular, it is square integrable. Since ipa is 
an eigenfunction of a symmetric difference operator whose coefficients have linear 
growth at ±oo, and since to different indices a correspond different eigenvalues, we 
conclude that these functions are pairwise orthogonal in i'^{Z'). 

Let us show that ||V'o|P = 1- Write (2.5), where we set x = y. Then the whole 
expression simplifies because (2.6) becomes equal to 1. Next, in the double contour 
integral, we replace the variable lu2 by its inverse. We obtain 

{M^;z,z',^)r = l-^^^' (l-^^^.y'^"-^ (l- VC'^r^)'"""' 
V / V / V ^^2 / wi a;2 
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To cvahiatc the squared norm we have to sum this expression over x € Z'. We spUt 
the sum into two parts according to the sphtting U = Z'_ U Z'_^_. We take as the 
contours concentric circles such that |wi| < \co2\ in the sum over Z'_, and |a;i| > \cu2\ 
in the sum over This gives us 

J J U;2-L0i LOi 102 J J LOl- ^2 ^2 

''^^ |wi|<|w2| |a;i|>|a;2| 

with 



(2TTi) 



Recall that both contours go in positive direction. 

Let us transform the second double-contour integral: keeping the second contour 
fixed we move tlic^ first contour inside the second contour. Then we obtain a double- 
contour integral which is cancelled out with the first double-contour integral, plus 
a single-contour integral arising from the residue of the function ui — > {u)i — U2)~^- 
l-^i' , .do; 1 - g / dw 



2m J ' ' ' oj^ 2-Ki J [l - ^Lo){ijj - ^) 

Thus, we have shown that the functions form an orthonormal family in i'^{U), 
and it remains to prove that this family is complete. For a; € Z', let 5x stand for the 
delta function at x. Since the functions 5x form an orthonormal basis, it suffices to 
check that 

{{5x,i}a)p{i.')f = ^{i}a{x;z,z\£,)f = 1, Vx G Z'. 

aeZ' aeZ' 

But this follows from the previous claim and the symmetry a a; established in 
the next proposition. □ 

Proposition 2.5. The following symmetry relation holds 
ipaix; z, z', ^) = i^xia; -z, -z', 0- 

Proof. Using the classical formula 



^ ^' cos(7r^) 

and the fact that both x + ^ and a + i are integers we check that 

r{z + x+^)T{z' + x+^) _ T{-z + a + ^)r{-z' + a+^) 

T{z-a+ i)r(y - a + i) ~ r(-^ -x+ ^)T{-z' - x + ^) ' 

Applying this to (2.1) and using another classical formula, 

F{A, B; C; w) = {1 - wf-^-^F{C -A,C-B;C; w), 

we get the required relation. 

Another way is to make a change of the variable in integral (2.4): 

CO CO = —= . 

This is an involutive transformation such that ^ ^-^d oo <-> l/VJ- As is 
readily verified, it leads to transformation (a, x, z, z') {x, a, —z, —z'). □ 
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Corollary 2.6. The functions tpa = i^aix; z, z' , ^) satisfy the following three-term 
relation 



(1 - Oxi>a = yj(,{z-a+\){z' -a+\)i)a-i 

+ ^S,{z-a-\){z' -a-\) xPa+i + {-a + £,{z + z' - aj) V-a- 

Proof. Under symmetry a; <-> a (Proposition 2.4), this turns into the formula stated 
in Proposition 2.1. Of course, a direct verification is also possible. □ 

The formulas of Proposition 2.1 and Corollary 2.6 show that the functions 
i)a{x; z, z' , possess the bispectrality property in the sense of [Gr]. 

Proposition 2.7. One more symmetry relation holds: 

Va(x; z, z', = (-l)"+>-a(-a;; -z, -z', 0, x,a€ U . 

Proof. This follows from the relation 

F{A, B; C; w) _ ^-^_c r{A-C + 1)T{B - C + 1) 



T{C) r(A)r(B) 

F{A -C+1,B-C + 1;2-C;w) 
r(2 - C) ' 



Cez, 



see [Er, 2.8 (19)]. Another way is to make a change of the variable, oj i— > l/w, in 
integral (2.4). □ 

In the remaining part of the section we will explain how the functions tjja are 
related to the Meixner polynomials. 

Let Z+ = {0,1,2,...}. Elements of Z+ will be denoted by symbols x,y (wc 
reserve the notation x and y for points of the lattice Z'). Recall that the Meixner 
polynomials are the orthogonal polynomials with respect to the weight function 

^! = r(a + l)^! ' 

on Z+, where a > —1 and, as before, ^ e (0,1). Our notation for these poly- 
nomials is ^fftm{x]a^,^). We use the same normalization of the polynomials as in 
the handbook [KS] (there is only a minor difference in notation: our parameter a 
corresponds to parameter (3 = a + 1 in [KS] , while our ^ is precisely parameter c 
in [KS]). 
Set 

mn{x;a,^) 



OT„(i;a,e) = (-l)"p^^^77^^1^a,«(i), ieZ+, (2.8) 
where 

■OC' 

The factor (—1)" is introduced for convenience: it will compensate the same factor 
in formula (2.10) below. 
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Proposition 2.8. Drop the assurnption that {z, z') is not in the degenerate series, 
and assume, just on the contrary, that z — N and z' = N + a, where N = 1,2,... 
and a > — 1. Then expression (2.1) for the functions ipaix; z, z', ^) still makes sense 
provided that 

x:=x + N-^, n:=N-a-\ (2.9) 
are in Z+, and in this notation we have 

ipaix; z, z', = 9n„(i; a, 0- 



Proof. As is well known, the Meixner polynomials can be expressed through the 
Gauss hypergeometric function in two different ways: 

9Jl„(x; a, = F{-n, -x; 

and 

r(x + i)r(-a-n) fi-^Y -^(-"'-"-";* + i-";|3t) 



mi„(i;a,0 = 



r(-a) V C y T{x + l-n) 

(-l)"r(a + l)r(a+ 1) fl-^y F{-n, -a-n;x + l-n; ^) 



T{a + n+l) \ ^ J T{x + l-n) 



, (2.10) 



see [KS]. Although the first expression looks simpler, it turns out that only the 
second expression is suitable for our purposes. Note that 

mi -a oil- = e(l-0"+^r(a + n+l) 
\WJ^n[ ,Q,UII r(a + l)r(n+l) 

Prom the last two formulas and the definition of Tin we obtain 



Vr(n + l)r(n + a + l)^ U 4 J 

F{—n, —a — n;x + I — n; ^^j) 
^ r(i + 1 - n) ■ 

Comparing this with (2.1) and taking into account (2.9) we get the required equal- 
ity. □ 

Thus, our functions 2pa can be obtained from the Meixner polynomials by the 
following procedure: 

• We replace the initial polynomials $H„ by the functions Tin ■ This step is quite 
clear: as a result we get functions which form an orthonormal basis in the space 
on Z+ with respect to the weight function 1. 

• Next, we make a change of the argument. Namely, we introduce an additional 
parameter N = 1, 2, . . . and we set a; = x— A^+^. Then we get orthogonal functions 
on the subset 

{-N + i, -AT + |, -TV + |, . . . } c Z' 
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which exhausts the whole U in the hmit as TV goes to infinity. 

• Then we also need a change of the index. Namely, instead of n we have to take 
a = N — n—\. We cannot give a conceptual explanation of this transformation, it 
is dictated by the formulas. Again, the range of the possible values for a becomes 
larger together with A^, and in the limit as ^ +00 we get the whole lattice U . 

• Finally, we make a (formal) analytic continuation in parameters N and a, 
using an appropriate analytic expression for the Meixner polynomials. Note that 
the difference equation of Proposition 2.1 and the three-term relation of Corollary 
2.6 precisely correspond to similar relations for the Meixner polynomials. 

We hope that this detailed explanation will help the reader to perceive the ana- 
lytic continuation arguments in Sections 3 and 7. 

Of course, instead of the lattice Z' we could equally well deal with the lattice 
Z, and then numerous "5" would disappear. However, dealing with the lattice 11 
makes main formulas more symmetric. 

3. The discrete hypergeometric kernel 

Let X be a countable set. By a point configuration in X we mean any subset 
X C X. Let Conf(X) be the set of all point configurations; this is a compact 
space. Assume we are given a probability measure on Conf(X) so that we can 
speak about the random point configuration in X. The nth correlation function of 
our probability measure (where n= 1,2,...) is defined by 

Pn{xi, . . . , Xn) = Prob{the random configuration contains xi, . . . , Xn}, 

where xi, . . . ,Xn are pairwise distinct points in X. The collection of all correlation 
functions determines the initial probability measure uniquely. 

We say that our probability measure is determinantal if there exists a function 
K{x,y) on X X X such that 

Pn{xi, . . . ,Xn) = det[K{xi,Xj)]'^.^^ , n = l,2,... (3.1) 

It is worth noting that if such a function K{x,y) exists, then it is not unique. 
Indeed, any "gauge transformation" of the form 

K{x,y)^^K{x,y), (3.2) 

f{y) 

where / is a nonvanishing function on X, does not affect the determinants in the 
right-hand side of (3.1). 

Any function K{x,y) satisfying (3.1) will be called a correlation kernel of the 
initial determinantal measure. Two kernels giving the same system of correlation 
functions will be called equivalent. 

As in §2, we arc dealing with the lattice Z' of (proper) half-integers. We split it 
into two parts, Z' = Z'_ U Z^, where Z'_ consists of all negative half-integers and 
1', consists of all positive half-integers. For an arbitrary A € Y we set 



X(A) = {Ai-i + i |i = l,2,...}cZ'. 
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For instance, ^(0) = The correspondence A ^(A) is a bijection between 
the Young diagrams A and those (infinite) subsets X dU for which the symmetric 
difference XA Z'_ is a finite set with equally many points in Z'_j_ and Z'_ . Note that 

X(A') = -(Z'\X(A)). 

Using the correspondence A i— > X (A) we can interpret any probability measure 
M on Y as a probability measure on Conf(Z'). This makes it possible to speak 
about the correlation functions of M. Our goal is to compute them explicitly for 
the z-measures. 

Now we can state the main results of the section.^ 

Theorem 3.1. For any admissible pair of parameters {z,z'), see §i, the corre- 
sponding mixed z-measure Mz,z',i * determinantal measure. 

Theorem 3.2. // {z,z') is not in the degenerate series (so that z and z' are not 
integers) then the correlation kernel of M^^^i^^ can he written in the form 

Kz,z',dx,y) = ■4>a{x;z,z',^)tpa{y;z,z',^), x,yez', (3.3) 

where the functions tjja are defined in (2.1). 

Note that the series in the right-hand side is absolutely convergent. Indeed, 
since {tJja} is an orthonormal basis in £^(Z') (Proposition 2.4), this follows from the 
fact that the series can be written as 

where stands for the delta-function at point x on the lattice Z', and (•, •) 
denotes the inner product in £^(Z'). 

Formula (3.3) simply means that K^z,z' ,£,{x,y) is the matrix of the orthogonal 
projection operator in i'^ilJ) whose range is the subspace spanned by the basis 
vectors V'a with index a e Z^ C Z'. 

Theorem 3.3. The correlation kernel (3.3) can also he written in the form 

Kz^z',dx,y) = 'Pz,z'{x,y)Kz^z>^^{x,y) (3.4) 

where, as in (2.6), 

Jr{x + Z + i)r(.T + z' + i)r(y + z+ ^)T{y + z' + ^) 
T{x + z' + ^)T{y + z + ^) 



^These results are essentially not new, see [B02], [BOS], [Ok2], [BOk, Example 3], and the 
comments at the end of the section. However, the method of proof is new. The same method, with 
suitable modifications, is applied in §7 for the computation of the dynamical correlation functions. 
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and 

1-e 



(27ri)2 

{wi} {uJ2} 



(i-^..r^(i-|)'(i-v?..)-(i-|)' 

J J wia;2 - 1 



1 1 

-2, ry-2 



X rfwi da;2 (3.6) 



where {coi} and {1^2} are arbitrary simple contours satisfying the following three 

conditions: 

• both contours go around in positive direction; 

• the point ^^/^ is in the interior of each of the contours while the point 
lies outside them; 

• the contour {wf^} is contained in the interior of the contour {u)2} {equivalently, 

{w2^} is contained in the interior of {uji}). 

The kernels K f (x,y) and I^z,z,i{xjy) O'l^e equivalent. Namely, they are re- 
lated by a "gauge transformation", 

Kz,z',dx,y) = -(^'^'^^1 Kz,z',dx,y), x,y e U, 
jz,z'\y) 

where 

^nx + z+\)V(x + z' + \) 

The kernel K_z,z,i{x,y) can serve as a correlation kernel for all admissible values 
of parameters {z,z'), including the degenerate series. 

Proof of Theorems 3.1-3.3. We prove these three theorems simultaneously. Let 
p{z,z _ _ _ ^2;„) denote the n-point correlation function of Mz,z',i- The proof 

splits into two parts. 

In the first part, we compute p^f '^^ for special values of the parameters (the 
degenerate series): 2; = = 1, 2, . . . and z' = N + a, where a > —1. Here we 
use the fact that for such {z,z'), the mixed z-measure can be interpreted as the 
so-called A''-particle Meixner ensemble. We show that formula 

pi''''^^\xi, . . . ,.x„) = det[K,^,,^dx^,XJ)]l^^, 

is valid (in particular, the values of the kernel in the right-hand size are well defined) 
when z = N, z' = z + a, provided that N is so large that the numbers Xi+ N — ^ 
are nonnegative. Then we check that in that formula, the kernel ^Kz,^',^ can be 
replaced by the kernel ^z,z',(- 

p(f'^''«Ha;i, ...,Xn) = det\Kz,z',dxi,Xj)]lj^i 

In the second part, we extend the latter formula to arbitrary admissible (z, z'). 
To do this we show that both sides are analytic functions in parameters {z,z',^). 
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Moreover, these functions are of such a kind that they are uniquely defined by their 
values at points {z = N, z' = N + a,^). 
We proceed to the detailed proof. 

Let, as in §2, iV be a natural number and a > —1. Consider the Meixner 
weight function Wa.^ on Z+, sec (2.7). The A^-point Meixner ensemble is formed 
by random A^-point configurations X — (ii >••• > xn) inX = Z_|_, where 

n 

Prob(X) = const • JJ {xi - xjf JJ Wa,^(ij). 

l<i<j<N i=l 

By the A^th Meixner measure we mean the corresponding probability measure on 
Conf(Z+). 

Lemma 3.4. The Nth Meixner measure is a determinantal measure. As its cor- 
relation kernel on Z+ x Z+ one can take the "Meixner kernel" 

N-l 

y)=J2 ^'»(^; ^m{y; a, 0, i,ye Z+, 

m=0 

where the functions £!Jtm(x; a, ^) are defined in (2.8). 

Proof. This is a special case of a well-known general claim about orthogonal poly- 
nomial ensembles, see, e.g., [De2]. □ 

Let Y(iV) C Y denote the set of Yoimg diagrams A with £{X) < N. Recall 
that the mixed z-measure with parameters z = N, z' > N — 1 is concentrated 
on Y{N). We define a bijection between Young diagrams A e Y{N) and A/'-point 
configurations X C Z+ as follows 

Xh^ X{X) = {xi, . . . ,xn), Xi = Xi-i + N, i = l,...,N. 

Lemma 3.5 ([B02, Proposition 4.1]). The correspondence X i— > X{X) takes the 
z-measure M^^z',^ with parameters z = N,z' = N + ato the Nth Meixner measure 

with parameters a, ^. 

Proof. Direct verification. □ 

Recall that we identify Mz^z',e. with its push-forward under the correspondence 
X^X{X). 

Corollary 3.6. Let z = N — 1,2, . . . and z' = z + a with a > —1. Assume that 
xi,. . .,Xn lie in the subset Z+ — AT + i C Z', 50 that the points Xi := Xi + N — ^ 
are in Z+ . 
Then 

pt^''i\xu ...,xn)= det [K^Zri^i' 

Proof. Let A € Y{N). Comparing the definition of the infinite configuration X (A) C 
Z' with that of the A''-point configuration X{X) we see that 

X{X) = (X(A)+Ar- i)nZ+. 

Then the claim follows from Lemmas 3.4 and 3.5. □ 

We take (3.3) as the definition of the kernel K_z,z',i{x,y). 
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Lemma 3.7. Let z ^ N = 1,2,... and z' = z + a with a > — 1. Assume that x 
and y lie in the subset — N + ^ C Z' , so that x := x + N— ^ and y := y + N — ^ 
are in Z+ . 

Then expression (3.3) for the kernel K_z,z',i{x,y) is well defined and we have 



Proof. We have to prove that 



^ ^pa{x;z,z',^)iljaiy;z,z',^) = ^ Ttmix; a, ^) Tlm{y; a, ^) (3.6 



m=0 



We recall that the functions tpaix; z, z', ^) were defined under the assumption that 
both z,z' are not integers. However, as it can be seen from (2.1), each summand 
in the left-hand side of (3.8) makes sense under the hypotheses of the lemma. 

Set 

a{m) = N -m - ^, m = 0,1, . . . ,N - 1 



By Proposition 2.8, 
which implies that 

JV-i _ 

J2 Mx;z,z',OMy;z,z',o = 'xiim{x;a,o^m{y;a,^) (3.9) 



m=0 



Finally, observe that 
1 



T{z-a+^) 



1 AT^i r(7v-a + i; 



a=Ar+i JV+I,., 



Wo conclude that the infinite sum in the left-hand side of (3.8) actually coincides 
with the finite sum in (3.9). □ 

Together with Corollary 3.6 this implies 

Corollary 3.8. Let z = N = 1,2,... and z' ~ z + a with a > —1. Assume that 
Xi, . . . ,Xn lie in the subset Z_|_ — + ^ C Z', so that the points Xi := Xi + N — ^ 
are in Z_|_ . 
Then 

Pn'^^'^Kxi,. ..,Xn)= det \I^z,z',d^i, Xj)]"^^^ 
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Lemma 3.9. Assume that 

• either {z, z') is not in the degenerate series and x,y €Z' are arbitrary 

• or z = N = 1,2, ... , z' > N — 1, and both x, y are in Z+ — N + ^. 

Then the kernel K_z,z' ,^{x,y) of Theorem 3.3 is related to the kernel K_z,z' ,^{x,y) 
by equality (3.4). Equivalently, the kernels are related by the "gauge transformation" 
(3.2), 

Kz',dx,y) = jf^K,,.-,dx.y\ (3.10) 

where fz,z' is defined in (3.7). 

Proof. Let us start with expression (3.3) of the kernel K_z.z'.(_ let us replace 
each sunimand by its integral representation (2.5). It is convenient to set a — ^ = k 
so that as a ranges over Z'_^_, k ranges over {0, 1,2,...}. Then we obtain 



(27ri) 



Kz,z',i{x,y) = (pz,z'{x,y) 

/ /(i-v/«™r'(i-f)'"(i-v«-.)-(i-^ 



2 Z-^ 



^ ^ ^(1 — v^a;i)(l — v^W2)\'^ dwi (iw2 



We can choose the contours {coi} and {0^2 } so that they are contained in the domain 
> 1. Since the fractional-linear transformation 

preserves the unit circle |w| = 1 and maps its exterior > 1 into its interior 
|c<;| < 1, we have on the product of the contours a bound of the form 



(wi - \/f)(a;2 - V^) 



<q<l. 



Therefore, we can interchange summation and integration and then sum the arising 
geometric progression in the integrand: 



E 

fe=0 



(1 - V?a;i)(l - V?W2)\ V ^1 



k 1 1 W1W2 



UJl 



(wi- V^)(a;2-v^) ; (1 - 0(c^i'^2 - 1) 



Then we obtain equality (3.4) with integral (3.6), as desired. Finally, wc can relax 
the assumption on the contour: it suffices to assume that {w^^} is strictly contained 
inside {1^2}^ as in the formulation of Theorem 3.3. 

It rcuirtius to show that (3.4) is equivalent to (3.10). According to (3.5) consider 
the expression 

1 r{x + z'+'.)Tiy + z+'.) 



^z,z' {x,y) 



^V{x + z + \)T{x + z' + \)T{y + z + \)T{y + z' + \) 
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Let US show that 

'Pz,z'{x,y) fz,z'{y) 

Indeed, l/^Pz,z' has the form 

a{x)b{y) 

and our hypotheses imply that a{x)b{x) and a{y)b{y) are real and strictly positive. 
We also have 

f ,(:,) "("0 _ 

JZ,Z'\X) — y ■ 

^ya{x)b{x) 

Therefore, we get 

fz,z'{x) ^ a{x)^Ja{y)b{y) ^ a{x)a{y)h{y) 
fz,z'{y) y/aix)b(x) a{y) a{x]b{x)a{y)h{y) a{y) 
^ 0'{x)b{y) ^ 1 
^Ja{x)b{x)a{y)b{y) ipz,z'{x,y) 

□ 

Corollary 3.10. Let z = N = 1,2, .. . and z' > N -1. Then 

n 

(3.11) 

provided that all the points xi, . . . ,Xn & lie in the subset Z+ — N + ^ c Z'. 

Proof. Indeed, this follows from Lemma 3.9 and Corollary 3.8. □ 

This completes the first part of the proof. Now we proceed to the second part. 

Lemma 3.11. (i) Fix an arbitrary set of Young diagrams V CY. For any fixed 
admissible pair of parameters {z,z'), the function 

xeT> 

which is initially defined on the interval (0, 1), can be extended to a holomorphic 

function in the unit disk |^| < 1. 

(ii) Consider the Taylor expansion of this function at = 0, 

oo 

xeT> fc=o 

Then the coefficients Gk,v{z, z') are polynomial functions in z, z'. That is, they are 
restrictions of polynomial functions to the set of admissible values {z,z'). 



Jz,z',i) 



(a;i,...,a;„) = det K^^^,^^{xi,Xj) 
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Proof, (i) Set D„ = P n Y„. By the definition of M^,^/,^, 



n! 

n=0 \AeX)„ / 

Eacii interior sum is nonnegative and does not exceed 1. On tiie otiier iiand, 

Ti=0 n=0 

This proves the first claim, 
(ii) By (1.11), 



E M,,,,,5(A) = (1 - e)^^' 13 13 (^)A(/)Af 



( dim A ^ ^ 



V n! 



It follows that 

n=0 ^ Ael>„ ^ ' 

Since each Vn is a finite set, this expression is a polynomial in z, z'. □ 

Now we can complete the proof of the theorems. Fix n and an arbitrary n point 
subset X = {xi, . . . ,x„} C Z', and regard pn'^ '^\xi, . . . ,x„) as a function of 
parameters z, z', ^. We want to show that equality (3.11) holds for any admissible 
{z, z'). Apply Lemma 3.11 to the set V of those diagrams A for which X (A) contains 
X, and observe that 

p(f'^''«)(a;i,...,ar„) = EM,,,,,5(A). 

xeT> 

It follows that pn'^ ■ ■ ■ ; Xn) is a real-analytic function of ^ S (0, 1) which 

admits a holomorphic extension to the open unit disk |^| < 1. Moreover, the 
Taylor coefficients of this function depend on z,z' polynomially. 

On the other hand, from the expression (3.6) for the kernel if 2. ^/^^ (a;, y) it follows 
that this kernel (and hence the right-hand side of (3.11)) has the same property, 
with ^ replaced by 

Thus, both sides of (3.11) can be viewed as (restrictions of) holomorphic func- 
tions in with polynomial Taylor coefficients. Since the set 

{{z,z') I z is a large natural number A'' and 2:' > A' — 1} 

is a set of uniqueness for polynomials in two variables, we conclude that equality 
(3.11) is true for any admissible {z,z'). 
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This proves Theorem 3.1 and Theorem 3.3. Now, Theorem 3.2 follows from 

Theorem 3.3 and Lemma 3.9. □ 

Comments. 1. The correlation functions of the z-measures M^^z',^ were first com- 
puted in [B02] in a different form: in that paper we dealt with another embedding 
of partitions into the set of lattice point configurations (in the notation of §8, we 
used the map X<-* X{X), instead oi X<—^ 2L (A)). The kernel Kz,z',i{3;, y) coincides 
with one of the "blocks" of the kernel considered in [B02]. The relation between 
both kernels is discussed in detail in [B05] (see also §8 below) . The proofs in [B02] 
and [BOS] are very different from the arguments of the present section. 

2. Two other derivations of the kernel K_z^z',^{x,y) are given in Okounkov's 
papers [Ok2] and [Okl]. In both these papers, the correlation functions are ex- 
pressed through the vacuum state expectations of certain operators in the infinite 
wedge Fock space. A (substantial) difference between the methods of [Ok2] and 
[Okl] consists in the concrete choice of operators. The general formalism of Schur 
measures presented in [Okl] is complemented by explicit computations in [BOk, 
§4]. 

3. As shown in the papers listed above, the kernel K_z,z'4{x,y) can be written 
in the form 

X — y 

where P and Q arc certain functions on Z' depending on parameters z, z', ^. Since 
P and Q are expressed through the Gauss hypergeometric function, we called 
K_z,z' ,(,{x,y) the discrete hypergeometric kernel. In general, kernels admitting such 
an expression are called integrable kernels, in accordance with the terminology of 
[IIKS], [Del], [B2]. 

4. The integrable form (3.12) can be readily derived from (3.3) using the three- 
term relation for functions tjja{x) = il>a{x; z, z'^) given in Corollary 2.6. Specifically, 
we obtain 

Kz,z',,{x,y) = ^ l-y ' (3.13) 

This derivation of (3.13) from (3.3) is quite similar to the standard derivation of 
the Christoffel-Darboux formula for an arbitrary system of orthogonal polynomi- 
als. Since, as explained in §2, the functions Va are closely related to the Meixner 
polynomials, this analogy is not surprising. 

5. Once we know that the functions form an orthonormal basis (Proposition 
2.4), the series expression (3.3) for the Vevnel K_z^z',(,{x, y) immediately implies that 
it is a projection kernel. This fact was first proved in [B05, §5] in a different way. 

6. The series representation (3.3) is equivalent to formula (3.16) in [Ok2]. A 
double contour integral representation of various correlation kernels related to Schur 
measures appeared earlier in [BOk] . 

4. Construction of Markov processes 

The goal of this section is to explain the construction of the continuous time 
Markov processes on partitions which will be studied in the rest of the paper. Their 
fixed time distributions are the z-measures considered in the previous sections. 

It is fairly easy to give the jump rates for these processes. However, it is not a 
priori clear why these rates define the process uniquely. Since we were unable to 
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find suitable uniqueness theorems in the literature, we will actually prove that the 
rates define the process uniquely and compute the transition probabilities for an 
underlying birth-death process. 

4.1. Preliminaries on Markov processes. Let us recall some basic facts about 
continuous time Markov processes and introduce the notation. 

The time parameter t always ranges over an open interval {tmimtmax) where 
tmin e K U {—00} and tmax € R U {+00}. Let us denote the state space by A, it is 
assumed to be either finite or countable. 

We also denote by P{s,t), s < t, the matrix of transition probabilities of a 
Markov process. This is a matrix with rows and columns marked by elements 
of A, its elements will be denoted by Pab{s,t), a,b E A. By definition, Pab{s,t) 
is the probability that the process will be in the state b at the time moment t 
conditioned that it is in the state a at time s. Thus, all matrix elements of P(s, t) 
are nonncgativc, and its sum is equal to one along any row. Such matrices are called 
stochastic. The transition matrices P{s,t) also satisfy the Chapman- Kolmogorov 
equation 

P{s, t)P{t, u) = P{s, u), s<t<u. (4.1) 

We assume that there exist A x A matrices Q{t) with continuously depending 
on t entries, such that 

Pab{s, t) = 5ab + Qab{t){t - s) + o{\t - s\), \t - s\ ^ 0. 

This relation implies that Qab{t) > for a ^ 6 and Qaa{t) < 0. Further, we assume 
that 

J2 Q'^bit) = -Qaa{t), for any a G A. (4.2) 

bj^a 

This is the infinitesimal analog of the condition J2beA-Pab{s,t) = 1. 

It is well known that (4.1) then implies that P{s,t) then satisfies Kolmogorov's 
backward equation 

-^P(s,t)=Q(s)P(s,i), s<t, (4.3) 
with the initial condition 

P{t,t)= Id. (4.4) 

Under certain additional conditions, P{s, t) will also satisfy Kolmogorov's forward 
equation 

^^P{s,t)=P{s,t)Q{t). (4.5) 

In our concrete situation we would like to define a Markov process by specifying 
the transition rates Q{t). However, it may happen that this does not specify the 
process uniquely (then the backward equation has many solutions P(s, t)). Unique- 
ness always holds if A is finite or, more generally, if A is infinite but the functions 
|Qaa(i)| are bounded on any closed time interval (see, e.g. [Fel]). However, these 
conditions are not satisfied in our case. There exist other, more involved uniqueness 
conditions for time homogeneous (stationary) Markov processes. However, in our 
approach, even if we restrict our attention to stationary processes, we still need to 
handle some non stationary processes as auxiliary objects. For these reasons we 
had to find some more special uniqueness condition. 
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Let us write Q{t) in the form Q{t) = ~R{t) + Q{t), where —R{t) is the diagonal 
part of Q{t) and Q{t) is the off-diagonal part of Q{t). In other words, 

nab\tj = —OabWaa[t), WabV') = i „ 

to, a = 0. 

For s < t set 

F{s, t) = exp R{T)dTj , G{s, t) = F{s, t)Q{s, t). 

Define Pl"l(s,t) recursively by 

p[°l(s,i) =F(s,i), P^''\s,t)= f G'(s,T)p["-il(r,i)dT, n>l, 

J s 

and set 

oo 

P(s,t) = ^pN(s,t), s<t. 

n=0 

Theorem 4.1 [Fel]. (i) The matrix P{s,t) is substochastic {i.e., its elements are 
nonnegative and J2b^>ib{s,t) < 1). Its elements are absolutely continuous and 
almost everywhere differentiate with respect to both s and t, and it provides a 
solution of Kolmogorov's backward and forward equations (4.3), (4.5) with the initial 
condition (4.4). 

(ii) P{s,t) also satisfies the Chapman- Kolmogorov equation (4.1). 

(Hi) P{s,t) is the minimal solution o/ (4.3) (or (4.5)) in the sense that for any 
other solution P{s,t) of {4.3) {or (4.5)) with the initial condition (4.4) in the class 
of substochastic matrices, one has Pab{s, t) > Pab{s, t) for any a,b G A. 

Corollary 4.2. // the minimal solution P(s,t) is stochastic {the sums of matrix 
elements along the rows are all equal to 1) then it is the unique solution of (4.3) 
(or (4.5)) with the initial condition (4.4) in the class of substochastic matrices. 

Let us note that the construction of P{s,t) is very natural: the summands 
P^^'(s,t) are the probabilities to go from a to 6 in n jumps. The condition of 
P(s, t) being stochastic exactly means that we cannot make infinitely many jumps 
in a finite amount of time. 

Our next goal is to provide a convenient suflacient condition for P{s,t) to be 
stochastic. 

Fix s e {tmin,tmax) and G G A. Foi any finite X , X c A, a G X , we denote 
by Tg^ a, X the time of the first exit from X under the condition that the process is 
in a at time s. Formally, we can modify A and Q{t) by contracting all the states 
6 e A \ X into one absorbing state b with Q~j^^ = Q for any c G X U {b}. We 

obtain a process with a finite number of states for which the solution P{s,t) of 
the backward equation is unique. Then Tg^ a,x is a random variable with values in 
{s, +oo] defined by 

Proh{T,,a,x <t} = P^l{s,t). 
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Proposition 4.3. Assume that for any a € A and any s < t, e > 0, there exists a 
finite set X{e) C A such that 

Prob{T,,„,x(e) <t} <£• 

Then the minimal solution P{s,t) provided by Theorem 4-1 is stochastic. 

Proof. Consider the modified process on tlie finite state space ^(e) U {b} described 
above. Since its transition matrix P{s, t) is stochastic, 

J2 Pabis,t) = l-P^-f,{s,t)>l-e. 

bex(£) 

The construction of the minimal sohition as the sum of p["l's, sec above, immedi- 
ately imphcs that Pab(s. > Pab{s,t). Thus. '}2b -^ab{s,t) > 1 — £ for any £ > 0. □ 

4.2. An application to birth-death processes. A birth-death process is a 
continuous time Markov process on A = Z+ = {0,1,2,...} such that the rates 
Qmn{t) vanish if \n — m\ > 1. In other words, the process can make jumps only of 
size 1. Our assumption (4.2) means that 

~Qnn{t) = Rnnit) = Qn,n+lit) + Qn,n-l{t)- 

Proposition 4.4. Assume that for any closed segment [t',t"] C (tmimtmax) there 
exists a sequence {7ri}^o of positive real numbers such that i?„„(t) < 7„ for any 
t € and X^nTrT"'^ — minimal solution P{s,t) is stochastic. 

Proof. We will apply Proposition 4.3. Let us fix a G A = Z+. As A = A(£) we will 
take a set of the form {0, 1, . . . , n — 1} for a suitable n. Then Ts,a,x is the moment 
of the first arrival at n given that we start at a at the time moment s. To simplify 
the notation, set T„ = Tg,a,{o,i,...,n-i}- 

In order to estimate r„ we will compare our inhomogeneous birth-death process 
to the pure birth homogeneous process with transition rates Qn,n+i = 7n, Qn+i,n = 
0, for all n € Let r„ denote the time of reaching n given that we start at a 
at time s. Note that since this is a pure birth process, once the process leaves 
{0, 1, . . . , n — 1} it never comes back. 

It is known that for the pure birth process with rates 7„ the minimal solution 
is stochastic if and only if J^n^n^ = [Fel], [Fe3, ch. XIV, §8]. By our 

hypothesis, '^^In^ = ^-^^ may denote by P{s,t) the unique stochastic 
solution of the backward and forward equations. Clearly, 

oo 

PT0h{fn <t} = Y,Pab{s,t), 
b=n 

which tends to zero as n — > oo. Thus, the statement of this proposition will follow 
from Proposition 4.3 if we show that Prob{T„ < t} < Prob{T„ < t} for any t > s. 
For any n = a + 1, a + 2, . . . , set 



Fr,{t) = Prob{T„ < t}, K{t) = Prob{T„ < t}. 
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Wc will prove that Fn{t) < Fn{t) for all n > a using induction on n. 
Let us start with n = o + 1. Clearly, 

i?„+i(t) = l-e^^»(*-«) 

because the time of the jump a ^ a + 1 is exponentially distributed with parameter 
7a. On the other hand, the probability that our birth-death process will jump from 
a to either a — 1 or a + 1 before time t equals 

1 _ !l{Qa,a + l(t)+Qa.a-l{l))dl Raa{t)dt ^ 

Thus, Fa+i{t) < 1 - exp(- Raa{t)dt), and since Raa{t) < 7a by hypothesis, the 
estimate follows. 

In order to prove the induction step, note that for n > a + 2 we have 
F„{t)= f rfF„_i(T)Prob{T„-T„_i<i-r|T„_i=r}, 

J S 

Fn{t) = ^ rfF„_i(r) Prob{f„ - f„_i < i - T I f„_i = r}. 
Arguing exactly as in the case n = a + 1 above, we see that 



Prob{T„ - r„_i < t - r I T„_i - r} = 1 



e 



Prob{T„ - Tn-i <t-T\ Tn-i = r} < 1 - e-T"-i(*-^). 
Hence, it suffices to verify that 



f dFn-^{T)(l - e-^"-(*--)) < f dK- 

J s J s 



L(T)(l-e-^"-(*-^)). 



When we integrate by parts both sides of this inequality, we notice that the non- 
integral terms vanish (because Fn-\{s) = F„_i(s) = 0). Thus, we obtain the 
equivalent inequality 

7n-i^V„_i(T)e-^''-(*-^'dt < 7n-i^'i?„-i(T)e-^"-(*-^)rff 

which immediately follows from _F„_i(r) < _F'„_i(r). □ 

Remark 4.5. It is very plausible that Proposition 4.4 holds under the weaker 
assumption Qn,n+i(i) < In- However, the proof of such a statement would require 
additional considerations. 

4.3. Birth-death process associated with Meixner polynomials. From now 

on we restrict our attention to birth-death processes with 

Qn,n+i = a{t){c + n), Q„,„-i = Pit) n, (4.6) 

where a{t) > 0, (3{t) > are continuous functions on {tmin,tmax) and c > is a 
constant. Proposition 4.4 implies that for any process of this kind there exists a 
unique stochastic solution P(s, t) of Kolmogorov's backward and forward equations. 
By breaking the time interval into finitely many subintervals, we may as well assume 
that a{t) and f3{t) are piecewise continuous functions with finitely many points of 
discontinuity at which they have finite left and right limits. 

The negative binomial distribution Wc,^ on Z+ with parameters c > and ^ G 
(0, 1) is defined by 

7re,4(n) = (l-Cr%r, n = 0,l,2,.... 



n 



It will be convenient to interpret 7rc,{ as an infinite row- vector. 
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Proposition 4.5. Let ^{t) be a continuous, piecewise continuously differentiable 
function in t with values in (0, 1). Assume that £^{t) solves the differential equation 

t G {tmint^max)' (^*'^) 



Then the row vector T^c^it) solves Trc,^(s)P(-s,t) = TTc,£,{t) for any s <t. 

Proof. Let us differentiate 7rc_j(s)P(s, with respect to s and use Kolmogorov's 
backward equation. Collecting the coefficients of Pxy{s,t) in the yth coordinate 
gives 

/ c£ x£ 

TTc,iis){x)Pa=v{s, t) I " ^ "f ^ ^"^'^ ^ '''' 

I , ,, 7^c,g(s)(a:- 1) ^, 7rc^(^)(a^+l)\ 
- a(c + a; - 1) , . fi[x + 1) ^ ' . 

Simplifications show that this expression is zero for all m, n if (4.7) holds. The 
initial condition 7rc,|(s)P(s,t)|^^j= T^c,i(€) is obviously satisfied. □ 

Once wc have a family of distributions T^c.at) satisfying TTc.^{s)P{s,t) = 7rc.j(4), 
we can define a birth-death process by the matrix of transition probabilities P{s, t) 
(which is uniquely determined by the jump rates) and one-dimensional distributions 
7rc,{(t)- 

It is not a priori clear what is a convenient way to parametrize these processes. 
In particular, multiplying both a{t) and (3{t) by the same function of t leads only 
to a reparamctrization of time in our process. In order to eliminate this freedom, 
we will always use one specific choice of time in our processes which we call interior 
or canonical time of the corresponding process. The convenience of this choice will 
soon become clear. 

The interior time is uniquely determined by the condition that a{t) and P{t) are 
expressed through ^{t) by 

Evidently, these formulas imply (4.7). Moreover, for any {a{t), P{t),^{t)) satisfying 
(4.7), if a{t) and l3{t) do not vanish simultaneously, we can choose a new time 
variable T{t) with 

so that [f~^a{t{T)),f~^P{t{T)),^{t{T)) satisfy both (4.7) and (4.8) as functions in 

T. 

Thus, from now on we will parametrize our processes by continuous, piece- 
wise continuously differentiable functions ^(t) taking values in (0,1) such that 
— 2 (this condition is necessary to guarantee the nonnegativity of a 
and /?). Such curves £,(t) will be called admissible. Then the corresponding birth- 
death process is determined by jump rates given by (4.6), (4.8) and one-dimensional 
distributions TTc,^{t)- We will denote this process by A^c,j( )- 
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In other words, if we set A{t) — — i In^(t) then A{t) has to satisfy three condi- 
tions: A{t) > 0, for aU t; \A{t)\ < 1 for all t; and A{t) is continuous and piecewise 
continuously differentiable. 

In terms of A{t) it is convenient to single out important special cases: A{t) = 
const corresponds to the homogeneous birth-death process; A{t) = t + const corre- 
sponds to pure death processes; and A(t) = —t + const corresponds to pure birth 
processes. 

Note that in case of a pure birth process A{t) will hit zero in finite time which 
means that in terms of the canonical time parametrization, the process reaches 
infinity in a finite amount of time. 

The connection of the processes -?Vc,5( ) with Meixner polynomials discussed in 
the previous section is already obvious from the fact that the distributions tTc.^ 
coincide, up to a constant factor, with the weight functions Wc-i,j, see (2.7). Our 
next goal is to express P{s, t) in terms of the Meixner polynomials. We will use the 
notation (2.8). 

Proposition 4.6. The matrix P{s, t) of transition probabilities for the birth-death 
process A^c,|(-) has the form 

P.v{s,t) = ' Ve"(«-*)mi„(a;;c,e(s))OT„(y;c,^(i)), (4.9) 



c,«(s) 



{x) 



n=0 



where s <t and x,y 

Comments. 1. In the stationary case ^{t) = const this formula was derived by 
Karlin and McGregor [KMG2] as a part of a much more general formalism, see also 
[KMGl]. 

2. The formula implies that P{s,t) depends on the initial value ^(s), final value 
^{t) and the length t — s oi the time interval. However, P{s,t) does not depend on 
the behavior of the curve ^(•) inside this time interval, as one might expect. 

3. The simplicity of the factor e"^*~*^ is a consequence of our choice of the 
interior time of the process. 

4. Since 9Jto(a;;c, ^) = (7rc,j(a;)) 2 , the prefactor may be rewritten as 

^c,m(y)Y _ mi{y;c,^{t)) 



5. The formula implies that the kernel on Z_|_ x Z_|_ 



{x,y) ^ ^g"$m„(a;;c,C)an„(y;c,r?) (4.10) 



71=0 



takes nonnegative values for C, € (0, 1) and < q < min | yj^, y'^j- The bound 
on q follows from the inequality |^/^| < 2. 

Our proof of Proposition 4.6 consists of few steps. Let us denote the right-hand 
side of (4.9) by P,j,(s,i). 
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First, we show that P{s,t) satisfies Kolmogorov's backward equation. Since we 
know that there exists only one stochastic solution, it remains to prove that P(s, t) 
is a stochastic matrix. 

The fact that the sum of the matrix elements along any row is equal to 1 is 
obvious (only n = term gives a nonzero contribution due to orthogonality of 
nonconstant Meixner polynomials to constants). The fact that P{s,t) is always 
nonnegative is not so obvious. In order to prove that we explicitly evaluate P{s, t) 
in the cases of pure birth and pure death processes, and then show that in the 
general case, P{s,t) is always a product of a "pure death" and a "pure birth" 
transition matrices. 

Lemma 4.7. The following relations hold 

^/^ix+l){x + c) Mnix + 1; c, + ^^x{x + c- l)m„{x - 1; c, 

-{x{l + + cO^„{x; c, = -n(l - 0^n{x; c, 0, 

(4.11) 

2$ ( 1 - C) ^ (x; c, C) + vW+W^ OT„ (x + 1 ; c, C) 

- ^/ix{x + c - 1) SW„(a; - 1; c, ^) = 0. 

(4.12) 



Proof. Straightforward computation using 

d TLX 

^2_ajl„(a;;c,0 = — S!7l„_i(a; - l,c+ 1,^ 

oE, c 

and [KS, 1.9.5, 1.9.6, 1.9.8]. □ 

Proof of Proposition 4-6. First of all, we need to verify that P(s, t) (the right-hand 
side of (4.9)) satisfies the backward equation. This is the equality 

--g^Pxy{s,t) = Qxx{s)Pxy{s,t) + Qx,x+l{s)Px+l,y{s,t) + Qx,x-l{s)Px-l,v{s,t) 

(4.13) 

with Qxxis) = -Qx,x+i{s) - Qx.x-i{s) and 

The computation proceeds as follows. One substitutes the sum in the right- 
hand side of (4.9) into the needed equality (4.13) and collects the coefficients of 
3n„(?/; c, ^(t)) using the relation d/ds = ^{s) d / d{^{s)) . Each such coefficient has 
two parts: one of them docs not involve ^(,s) while the other one is equal to ^(s) 
times an expression not involving ^(s). It turns out that each of these parts vanishes, 
for the ffist part the needed relation is (4.11), and for the second part one uses 
(4.12). The details arc tedious but straightforward, and we omit them. 

As was mentioned before, it remains to prove that Pxy{s, t) is always nonnegative. 

Let us use the notation P^ {s,t), P^{s,t) for P{s,t) when we consider a pure birth 
or a pure death process (that is, ^(r) = 6^^°"^*+^^ or ^(r) = e'^™^*"^'^, respectively). 
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Lemma 4.8. (s, t) is the unique solutions of the backward equation for the pure 
death process with the initial condition P{t,t) = Id, and P^(s,f) is the unique 
solution of the forward equation for the pure birth process with the same initial 
condition. Furthermore, with the notation ( = ^(s), r] = ^{t), we have 




a;! 



x>y, 



{x-yy.y\ ' (4.14) 

x<y, 



y-x 



{C + X)y- 



x<y, 



{y-x)\ ' (4.15) 
x> y. 

Proof. Consider P^{s,t) first. Since Qx,x+i = 0, Kolmogorov's backward equation 
takes the form 

= T3f(^ {-PU',t)+PUyis,t)) ■ (4.16) 

If wc fix y then these differential equations can be solved recursively: we subse- 
quently find Po^y, Pi,yj -^2,2/) • • • ) using the initial conditions P^y{t,t) = S^y This 
shows that the backward equation for the pure death process has a unique solution. 
A straightforward calculation shows that the expression in the right-hand side of 
(4.14) satisfies this equation (with ^(s) = e<=°"^*-2«). 

The case of the pure birth process is completely analogous. □ 

Since for the pure death process the backward equation has a unique solution, 
we have just shown that P{s,t) — P-^{s,t) is the corresponding transition matrix. 

In order to make a similar conclusion for the pure birth process, we need to know 
that P{s, t) satisfies the forward equation. This fact can be proved directly using 
Lemma 4.7. It can also be reduced to the case of the backward equation as follows. 

Note that T:^^^(^s){x)Pxy{s,t) remains invariant under the changes 

s^-t, t^ -s, CM ^ f(r) := C(-t), x^y. (4.17) 
Thus, instead of computing §iPxy{s,t) we may compute 



U — ~t.V — ~'S 



with the new ^(t) obtained from ^(t) by the time inversion. Since we already know 
that Pxy{s,t) solves the backward equation, for ^(r) = e'=°"'**+^^ we obtain,^ cf. 
(4.16), 

r ^y 



2y ^c,au)(y) 



+ 



(-Pyxiu,v) + Py-i^xiu,v)y 



®The cirgument goes through for any admissible curve ^(•), it just becomes more tedious. 



MARKOV PROCESSES ON PARTITIONS 



37 



Using 7r^_|^(„) (y) = tt^ ^^^^ (j/- 1) • (c + y- and substituting u = -t, v = -s, 

we obtain the needed forward equation 

d_p 2{c + y)m p ... , 2{c + y-l)m p . . 

The conclusion is that in case of the pure birth process, P{s,t) satisfies the 
forward equation, and by Lemma 4.8 we have P{s,t) = P^ {s,t). 

The nonnegativity of P{s, t) for arbitrary admissible curves ^(•) follows from 

Lemma 4.9. Let ^(•) he an admissible curve and -^c,?( ) the corresponding birth- 
death process. Then for any s < t, P{s,t) is a product of P^{s,u) with ^(t) = 
g-2(r-s)+in£(s) P'^ {u,t) with ^(r) = e2(-^-*)+'"«*) for a certain choice of u. 
Specifically, u is determined from the continuity condition: 

g_2(„-s)+i„€(s) ^ g2(u-t)+in5(t) ^ ^ ^ t±l + In^(s) -In^ft) ^ ^^^^^^ 

Proof. This statement follows from the Chapman-Kolmogorov equation (4.1), which 
P(s, t) obviously satisfies due to the orthogonality of Meixner polynomials, and from 

the fact that P(s, t) does not depend on the specific form of the curve ^(•), sec Com- 
ment 2 after the statement of Proposition 4.6. Thus, we may just replace ^(t) by 
a continuous combination of e~^'^+'=°"^* and e^'^+'=°"^* and preserve ^(s), ^{t), and 
t — s. Note that the fact that u given by the formula above is between s and t 
follows from the inequality <2. □ 

Lemma 4.9 implies that Pxy{s,t) is always nonnegative, and this completes the 
proof of Proposition 4.6. □ 

Corollary 4.10. The process obtained from -A^c,{( ) by the time reversion is also of 
the form with ^(r) = ^{—t). 

Proof. -^c,?( ) is characterized by the fact that it is a Markov process with two- 
dimensional distributions 

Prob{iVc,«(s) = X, iVx,C(i) =y} = T^c^s) i^)Pxy{s,t). 

As was already mentioned above, the right-hand side of (4.9) multiplied by ^^^^{s) {x) 
is invariant with respect to (4.17). This implies the statement. □ 

Note that, in particular, time inversion turns our pure birth process into the pure 
death process and vice versa (essentially, we gave a proof of this fact before Lemma 
4.9), and the stationary process iVc,{ with ^ = const is reversible. This is well 
known; any stationary birth-death process with an invariant measure is reversible 
with respect to this measure). 

4.4. Markov processes on partitions. Our next goal is to extend birth death 
processes ■^c,5(-) to partitions in the following sense. We construct continuous time 
Markov processes on the state space Y (the set of all Young diagrams, see §1) 
parametrized by admissible pairs {z,z'), see §1, and admissible curves ^(•). The 
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projection of such a process on Z_|_ obtained by looking at the number of boxes of 
the random Young diagrams, coincides with A^^2',j(.). 

Let us fix a pair (z, z') of admissible parameters and set c = zz' > 0. Given an 
admissible curve ^(•), wc define the matrix Q of jump rates of our future Markov 
process ^z,z',i on Y by (set n = |A|) 



<9am(s) 



(c + n) 1 + 



2i{8) i-e(s) 



•J3^^,(n, A;n+ 1,/i), 



n 1 



-(c + 



2i{s) ] 1 - i{s) 



■p^{n,X;n- 1,/x), 



1 + 



1 



1 



M = A, 



2^(s)y i-c(s) V i-^(s)' 

(4.19) 

and Q\f_i = in all other cases. Here pl^, and are transition and cotransition 
probabilities from §1, see (1.9) and (1.1), and the expressions involving ^ come from 
(4.8). Note that under the projection Y — > Z_|_, A i-^ |A|, this matrix Q turns into 
the matrix of jump rates for iVc,^( )- 

Proposition 4.11. The minimal solution P{s,t) of Kolmogorov's backward equa- 
tion with the matrix Q defined above is stochastic. 

Proof. We apply Propositions 4.3, 4.4. In the proof of Proposition 4.4 it was shown 
that for any a G Z+ there exists a set of the form X = {0,1, . . . ,n — 1} such that 
the probability of exiting X during the time period from s to t with tlie initial state 
a is smaller than any given positive number s. This means that if we start at time 
s from A e Y with |A| = o then the probability of exiting Yq U Yi U • • • U Y„_i 
before time t is just the same as for the birth-death process and, hence, is less than 
e. Proposition 4.3 concludes the proof. □ 

Proposition 4.12 (cf. Proposition 4.5). For any s <t 



Z,z',i{t)-: 



(4.20) 



where P{s,t) is the transition matrix of Proposition 4.11, and M^^^i^^ is the mixed 
z-measure (1.11) viewed as a row-vector with coordinates marked by elements ofY. 

Proof. Since the formula obviously holds for s = t, it suffices to show that the 
derivative with respect to s of the left-hand side of (4.20) vanishes. Thus, it suffices 
to show that 

-jf M,,,,,5(,)(m) + J2 ^.,.',«(a)(A)QA^ = (4.21) 



for any fi gY. Recall that 

M^,^/,j(A) =M^"j(A)7rc,^(n) with c = zz', n = \X\. 

Substituting this relation into (4.21) we notice that we can perform the summation 
over A using (1.5) and (1.6). Factoring out M^^'^P{fx) leads to the formula which 
states that the derivative of 7rc,|(s)P(s, f) with P{s,t) being the transition matrix 
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for the birth-death process Nc^^, with respect to s vanishes. But this has already 
been proved in Proposition 4.5. □ 

We conclude that given an admissible pair {z, z') and an admissible curve ^(r), 
there exists a unique continuous time Markov process on Y with jump rates Q de- 
fined above and with one-dimensional distributions M^, j,/ j(t-). This Markov process 
will be denoted by A2^j,/^j(.). 

As for the birth-death processes, we single out three important special cases: 
the stationary process ^ = const, the ascending process ^(r) = e2T+const ^^^^ ^-^q 
descending process ^(r) = _ The projections of these processes on Z+ are 

the stationary birth-death process, the pure birth and the pure death processes, 
respectively. 

As in §1, for A e Y we denote by dim A the number of ascending paths in the 
Young graph leading from to A. More generally, we denote by dim(/x. A) the 
number of ascending paths in Y leading from /i to A; if there are no such paths we 
set dim(/Lt, A) = 0. Also, for /i, A G Y such that /x C A we set 

(i,j)G\\fj. 

where the product is taken over all boxes in A \ /i. 

Proposition 4.13 (cf. Lemma 4.8). The transition matrix of the descending 
process A^^zi^^^.) has the form 

\(i-vKJ W-Ov) {x-yV-yi dimA ' 



and the transition matrix of the ascending process A^, has the form 

pT -f)- f f V-C \ dim/i dim(A,^) 

^"^'^"Vl-C; U-Cy {y-xy.y\ dimA '^''^ 



(4.23) 

where C = C(s), V = ^{t), x = |A|, y = 

Proof. Let us consider the descending process first. It is immediate to check that 
the matrix J-*;[^(s, t) obtained from the transition matrix P^y{s, t) of the pure death 
process by 

Pi^(s,t)=P4(s,t) 

x^jj^a^) A; x-1, At^^-^-i))pHa;-l, At^''"^"^^ a;-2, /x(^-^-2)) . . .pi(y+l, yu^i^ y. At) 

where the sum is taken over all paths /i — ^i^^'' [i^^^ /i*^^"^-* ~ A from \i 

to A, satisfies the backward equation. All terms in the above sum are equal to 
dim /x/ dim A, and the number of terms is equal to dim(Ai, A). Together with (4.14) 
this implies (4.22). 

Similarly, for the ascending process one has 

X Y,v\.' (a;, A; a; + 1, \^^^)v\^, (a; + 1, \^^^ ; a; + 2, A^^)) • • • p\^, {y-\, \^y—^) ■ y, /,) 
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where the sum is taken over all paths A = A(°) / A^^^ / ■■■ / A(^-«) = /i from A 
to /i. Again, the product of transition probabilities does not depend on the path 
and it is equal to 

1 a::! dim^ 
(c + ?y! dimA ■ ^^^^ 

while the number of paths is equal to dim(A,/u). Together with (4.15) this gives 
(4.23). □ 

5. Transition matrix for integral values of z. 

Our main goal in this section is to obtain a formula for the transition matrix of 
the process in the case when ^; is a nonnegative integer. For z = \, the 

process A^^/^j.) coincides with the birth death process f^z' .c (because it lives on 
the Young diagrams with only one row), and our formula is reduced to (4.9). 

Fix i; = A/'e{l,2, ...}. In order for (z, z') to be an admissible pair, we must 
have z' G K and z' > N — 1. We will use the notation z' = N + a, a > —1. As 
before, we set c = zz' = N{N + a). 

As was mentioned in §1, the support of Mjv,Ar+a,{ consists of the Young diagrams 
with no more than N rows. It is convenient to parameterize such diagrams A by 
sequences of N strictly decreasing nonnegative integers {xi, . . . ,Xn), 

Xi = Xi+N-i, i = l,...,N. 

Given an admissible curve ^(•), set 

oo 

i;«,t(a;,y) =^e'=(^-*)OTfc(a;;a,e(s))^fe(2/;a,^(t)), x,yeZ+. (5.1) 
fe=o 

Theorem 5.1. Let X, n be Young diagrams with no more than N rows, and let 

{xi, . . . ,xn), {yi, ■ ■ ■ jJJn) be the corresponding sets of decreasing nonnegative in- 
tegers. For any admissible curve ^(•) the transition matrix of the Markov process 
^N,N+a,i{-) has the form 

P\i,{,s,t)=e 2 — det t;«,t(a;i,yj) (5.2) 

We will use the term Karlin-McGregor representation for this formula. 

Proof. The arguments follow the same pattern as in the proof of Proposition 4.6 
(which is a special case of this theorem) . The first step is to show that the right- 
hand side of (5.2) satisfies Kolmogorov's backward equation. After that we prove 
that this solution is stochastic. 
We will use the notation 

(r) 

x = {xi,...,xn), y = {yi,...,yN), £r = (0, • . • ,0, 1 ,0, . . . ,0), l<r<N, 

c = m, i = m, n = \x\ = f2x,- ^^^~^\ 

fk{-)=m{-;a, ^(s)), gk{-)=mk{-;a, ^{t)). 
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Also, denote the right hand side of (5.2) by Pxy{s,t). 
The formulas of § 1 imply 



M^,iv+a,e(t) (a*) V _ (1 - 77) § / r(y, + a + i)r(x, + 1) ' v{y) 



Mjv,iv+a,«(.)(A); <;^(i _ r{xj +a + l)r{y, + 1) J V{x) 

(5.3) 

where V{u) = ni<i<j<Ar(^j ~ '"j) stands for the Vandermonde determinant. 

Note that adding to A one box or removing from A one box is equivalent to 
adding to x or subtracting from x, where r is the row number of the box. We 
have, cf. (1.9), (1.1), 



PN,N+ai'"'^ x;n+l,x + er) = 



Xr + a + 1 V{X + Er) 



p^{n, x;n — l,x — Sr) = 



c + n yix) 

Xr V{x — Er) 

^ V{x) ■ 



The right-hand sides of these relations conveniently vanish exactly when xzLsr does 
not represent a Young diagram (two of the coordinates of a; ± are equal). Using 
(4.19) we can now write down the needed backward equation: 



"ds 



2e(s) / 1 - m 



2e(s)/ 1 -as) 



+ I 1 + ^^TTT I TTT 2^ Px+er,y{s,t) 



ns)j i-c(s)^ 



+ 1- 



iis) 



1 ^ 

1 - £(s) ^ 



Vix) 

XrV{x — Er) 



J 1 - e(s) fr{ Vix) 



Px — Sr,y{^J ^)- 



It is time to use the definition of Pxy{s,t). The Cauchy-Binet identity (see, e.g.- 
[Ga, ch. I, §2]) implies (in all determinants below the indices run from 1 to A'') 



det[vs,t{xi,yj)] = det 



J2^^'~'^'fki^i)9k{yj) 

fe=0 

e(«-*)(^^+-+'=-)det[/fe.(x,)]det[<?fe.(y,)]. (5.4) 



fcl>fe2>--->/siV>0 



Let us use this relation and (5.2) for Pxy{s, t), Px±er,y{s, t) in the backward equation 
and collect the coefficients of det[gki{yj)]- Factoring out 



M 



N,N+a 



M, 



1_ 

' g(^-t)(fci+-+fe^-^^^) 



N,N+a,e,{s) 
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we obtain (Vandermonde determinants cancel out quite conveniently) 

/ N{N -1) 



(fci + • • • + 



N 



det[/fc,(a;j- 



+ + ^(^^ + l)ixr + a+ 1)C det [fk, {{x + er)j 



2C ) 1- 

+ I ^ " ^ ) ^ VMxr + a)C det[/fc,((x - er)j)] . 



N 



(5.5) 

We claim that the part the left-hand side of this relation that does not involve 
derivatives ^(s) and — ^ equals the part of the right-hand side without i{s), and 
the part of the left-hand side with derivatives equals that of the right-hand side 
with ^(s). 

The part without derivatives gives 

C(c + n) + n + (1 - - (fci + • • • + fcjv)(l - det[/fc, {xj)] 

N 



= Z!(\/(^'- + l)(^'- + " + l)C det[/fe,((a; + e^)j)] 

+ y^Xr{xr + a)C det[/fc,((a; - er)j)]). (5.6) 



r=l 



The right-hand side of this equality can be rewritten as 

N 

J2 Z!^Sna/fe^(^,(a;i)---/fe^(„,(a;Ar) 

^ V {Xr + l){Xr + a + i)Cfa{r){Xr + 1) + \/ XrjXr + q)C fa{r){Xr - 1) 

fcr{r) {Xr) 

By (4.11) the last ratio equals Xr{l + C) + {a + 1)( — k„^r){^ — C); and the whole 
expression equals 

{{xi + ■■■+ xn){1 + C) + N{a + 1)C - (fci + • • • + kN){l - 0) det[/fc, {xj)], 

which is exactly the left-hand side of (5.6) with c = N{N + a) and n = ~ 
N{N - l)/2. 

The part of (5.5) with derivatives gives 

- §; det[/,,(:r,-)] = Y^i^i'^r + l)ixr + a + !)( det[/fe.((x + e.),-)] 

- \/ Xr{xr -\- a)C det[fki{{x - er)j)])- 
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The same operation with determinants and (4.12) show that this equahty holds. 

This concludes the proof of the fact that P\^(s,t) (the right-hand side of (5.2)) 
satisfies the backward equation for K^^z', £,{■)■ It remains to verify that the Y x Y 
matrix P\n{s,t) is stochastic. 

Let us check that Px^is, t) = 1 for any A S Y and s <t. 

Using elementary row operations on the Vandermonde matrix, we obtain 

V{u) = H {ui- Uj) = det[uf-^] = (-1)^^ det[9Jljv-i(%;a,0] (5.7) 

l<i<j<N 

jV(jV-I) 

where the sign (—1) 2 appears because the highest coefficient of dyik[u; a,^) is 
(-1)'=. Then (5.2) and (5.3) imply 

(iet[jN-i{Xj)\ 

Apply (5.4) and sum the result over all t/i > • • • > t/jv > 0. The Cauchy-Binet 
identity gives 



X] idki iVj )] det [qn-z iVj )] = det 

vi>--->yN>0 



^9ki {u)gN-j{u) 



.«=o 



(5.8) 



Orthogonality of Meixner polynomials means that Y^x) gk{'u)gi{u) = Ski- Hence, 
the last determinant equals 1 if fc^ = N — i for all i = 1, . . . ,N and vanishes 
otherwise. This gives the desired result. 

Finally, the nonnegativity of Pxij,{s,t) follows from Proposition 4.13 and 

Lemma 5.2 (cf. Lemma 4.9). For any admissible curve £,{■), the matrix P{s,t) 
given by the right-hand side o/(5.2) is the product of the transition matrix P^{s,u) 
for the descending process with ^(r) — e^^('^~**)+^"^('*) and the transition matrix 
P^{u,t) for the ascending process with ^(r) = e^('^~*)+'"^(*) with u given by (4.18); 

_s + t In^(s) - In^(t) 



The proof of this lemma is very similar to that of Lemma 4.9. The Chapman- 
Kolmogorov equation for P{s,t) is easily verified by means of (5.4), the orthogo- 
nality of Meixner polynomials, and the trick with the Cauchy-Binet formula used 
above. 

Since the transition matrices of the descending and ascending processes have 
nonnegative matrix elements (4.22), (4.23), the matrix elements of P{s,t) are also 
nonnegative, and the proof of Theorem 5.1 is complete. □ 

Corollary 5.3. The kernel (4.10) is totally positive. 

Indeed, Theorem 5.1 shows that the minors of this kernel are, up to positive 
factors, transition probabilities. 



44 



ALEXEI BORODIN AND GRIGORI OLSHANSKI 



§6. Analytic continuation 

The goal of this section is to employ the analytic continuation in the parameters 
z,z' in two directions. First, we show that the factorization of the transition matrix 
for into the product of transition matrices for descending and ascending 

processes is carried over from the case of integral z (Lemma 5.2 above) to the 
general case. Second, we prove an analog of Lemma 3.11 for finite-dimensional 
distributions of A^ This result will be used later to compute the dynamical 

correlation functions of Az^zi^^(^.y 

Proposition 6.1. The statement of Lemma 5.2 holds for arbitrary admissible pair 
of parameters {z,z'). That is, the transition matrix P{s,t) for A^_^/_j(.) is equal to 
the product P^{s,u)P^ {u,t) with the descending and ascending processes and the 
time moment u specified as in Lemma 5.2. 

Comments. 1. This statement together with Proposition 4.13 allows to write P(s, t) 
down rather explicitly. Namely, we obtain 

_ iC-eYiv-eni-r^r' /idimi. 

- ^i^^ _ nIdimA ^^'^^ 



(1 _ _ eyn^rn ^-^^^^ ^) ^^^^^^ ^) 

(/-m)!(n-m)! {z)^{z')^ ^ ' ' 



fie 

with the notation 

C = ^{s), v = ^it), 6' = e^-*V^, / = |A|, m=\n\, n=\iy\. 

The sum above ia actually finite; fj. ranges over all Young diagrams which are smaller 
than both A and f. 

2. The formula for P\f^{s,t) above after the multiplication by 



M,,,,,f (A) = (1 - Cr'C\z)x{z') 



dim^ 



becomes symmetric with respect to {X,Q < — > {v^rf). This shows that the time 
reversion of A2_j:/^j(.) is again a process of this form with new ^(r) = ^(— r), cf. 
Corollary 4.10. 

Proof. Set P{s,t) = P^ {s,u)P^ {u,t). Clearly, this is a stochastic matrix that 
satisfies the initial condition P{t, t) = Id. Thus, it suffices to verify that it satisfies 
Kolmogorov's backward equation 

-^Px,.{s,t) = Qxx{s)Pxx{s,t)+ J2 Qxnis)PnAs,t) + Qxn{s)P^,,.{s,t) 

with Q{s) given by (4.19). As we substitute the expression in the right-hand side 
of (6.1) for P{s,t), we see that both sides of the equality above, as function in z 
and z', have the form 



1 — T] 



zz 

X 



{a polynomial in z,z'}. 
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Since the equality has aheady been estabUshcd for (z, z') = {N, N + a) with N = 
1,2,... and a > —1, it must hold for arbitrary {z, z'). □ 

The statement that we prove next will be used in the derivation of the dynamical 
correlation functions later in the next section. 

Take an admissible pair of parameters {z, z') and an admissible curve ^(•), and 
consider the Markov process K^^^i^^f^.y Let t\ < t2 < ■ ■■ < t„ be arbitrary time 
moments. Set 

Proposition 6.1 (or formula (6.1)) and the fact that the process is Markovian 
imply that the finite-dimensional distributions 

Prob{A;j_;j/_j(.)(ti) =X{i); i = l,...,n} 

with given Young diagrams Ai, . . . , A„, depend on parameters ^i, . . . , ^„ and 7?i2, 
??23, • • • , ?7n-i,ra but they do not depend on the behavior of ^{t) inside the intervals 
{ti, ti+i). Thus, in order to compute these finite-dimensional distributions we may 
replace our process by a sequence of alternating descending and ascending processes: 
We start off at the time moment ti and go down till time 

ti+t2 , lna-ln6 
«i2 = ^- + ^ , 

the value of ^ at this moment is exactly 7712 . Then we go up till ^2 and then again 
down, etc. The time moments when we change directions are 

il < W12 < *2 < "23 < • • • < Un^n+l < tn 

with = {ti + ti+i)/2 + (In^j — ln^i+i)/4, and the values of ^ at these points 

are r]i^i+\. At tiS we switch from going up to going down, and at Wj^i+i's we switch 
from going down to going up. 

Fix arbitrary subsets Pi, 252, • • • , T^n] 2^12, ^^23, • • • , T^n-i,n of Y. Let us compute 
the probability that our new descending-ascending process hits all Vi, at 
the time moments ti, Wj,j+i, respectively. It is equal to the sum 

^--^'-Ci (^(1)) ^i(l), ^(1,2) "I2)^il,2), A(2) {U12M) ■ ■ ■ 

\ii)eVi,i=l,...,n 
/^(j>j+l)el'j,j+i> J = lv,n-1 

-^A(n-l),;i(n-l,Ti)(*"-l''""-l.")-'^/I(n-l,n), (^•^) 

Proposition 4.13 shows that this is in fact a function in ^i,r]i^i+i which we will 

denote by F{^,ri) with ^ = (^1, . . . ,^„), rj = (7712, . . . ,t/„-i,„). The parameters 
and rji^i+i may take arbitrary values between and 1 subject to the inequalities 

> ^12 < C2 > '723 < ^3 > • • • < ^n-l < Cn- (6.3) 
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Proposition 6.2 (cf. Lemma 3.11). (i) The function F{^,r]) is a real-analytic 
function in ^ and rj. 

(ii) The function e ^ F{e^;er]) which has been defined so far for e > 0, can be 
analytically continued to a nonempty disc of the form \e\ < const. The coefficients 
of the Taylor decomposition of this function at e = are polynomial functions in 
z, z'. 

Comment. This statement implies that the function F{£_,ri) viewed also as a func- 
tion in z, z', is uniquely determined by its values on the arguments rj, z, z') with 
arbitrary 5, 77, < ^j,r7jj+i < 1, satisfying (6.3) and (z, z') = {N, N + a) with 
N = 1,2, ... , and a > —1. This uniqueness will be used in the next section to 
extend certain formulas derived in the case of integral z to the general case. 

Proof In order to prove (i), by Weierstrass' uniform convergence theorem it suffices 
to check that the series (6.2) with M^^^/^j^ (A(l)) and all , P^ replaced by their 
expressions given by (1.11), (4.22), (4.23), converges absolutely and uniformly in 
^, rj varying in small discs around their values. It is more convenient to work with 
the case when all T>i and coincide with Y; clearly, the needed convergence of 

the restricted sum follows from that of the unrestricted sum. 

Set li = |A(?')|, rrii^i^i = |/z(i, i + 1)|. As seen from the proof of Proposition 4.13, 
the matrix elements P^^ and split into products of transition probabilities for 
pure death and pure birth processes of Lemma 4.8 and (co)transition probabilities 
on the Young graph. By (1.11), Mz^z',£,{^) is also a product of the negative binomial 
distribution 7rc,^(n) on Z+ and the probability distributions M^"', (A) on Y„'s. The 
probabilities related to the Young graph do not depend on ^ and rj. Thus, we can 
split the sum in (6.2) (remember that all and 'Dj^j+i's are equal to Y) into two: 
the outer sum is taken over all nonnegative integers Zi, . . . , /„ and mi2, . . . , m„_i,„ 
satisfying 

h > mi2 <h> m23 < • • • > m„_i,„ < Z„, (6.4) 

and the inner sum is taken over all Young diagrams \{i) e Y; . , /u(i, i + 1) e Y^ . . 
The relations (1.5), (1.6) applied to the z-measures imply that the inner sum is 
equal to 

7r..',«i(Zi)P/^,„^,(ti, wi2)Pi,„i,(wi2,i2) • • • 

• • • -P/„_i , m„_i,„ ""-l.")^m„_i,„, i„ in) (6.5) 

with P^ and given by Lemma 4.8. Thus, we need to verify the uniform conver- 
gence of the sum of such products taken over nonnegative integers satisfying (6.4). 

Denote 

i = 1, . . . ,n, with 770,1 := 0, 
i = 1, . . . ,n — 1. 

The inequalities (6.3) imply that Ui G (0, 1) and u^.i+i G (0, 1] for all i. Clearly, 
small (complex) variations of ^ and ry lead to small (complex) variations of Wj's and 
u,,i+i's. What is important for us here is that if the variations are small enough 
than Uj's are bounded away from 1 and Uj^j+i's are bounded away from 0. 



_ ii - 

1 - Vi-i,i 

^ (1 - ii)Vi,i+i 

'''^^ (l-??M+l)6 
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Take the absolute value of (6.5) and sum it over /„. We have 



E 



X E 



- mn-l,n)l„- 



= ll - w„ 



{In 



^n— l,n 



\ZZ +mn-l,Ti 



1 



(6.6) 



We conclude that this expression can be estimated, cis a. function of i.m 

by a 

constant times a geometric progression of tlie form r™"-!-" with a suitable r > 0. If 
the variation of ?/„ is small enough, it is close to the real axis. Hence, by decreasing 
the variations we can take r arbitrarily close to 1. 

Let us substitute this estimate into (6.5) and sum over m„_i,„. We obtain 



const 



E 



(^n— 1 ) '^n—l,n) 



= const \Vn-l,n\ 



E 



= const 



1 



1 



(^n— 1 ^n— l,n) *^7i— l,n* 



-l,n| 



'^n— 1,71 

Again, this is bounded by a constant times a geometric progression r^n-i where r 
can be made arbitrarily close to 1 by considering small enough variations of u„ and 

^71—1,71* 

The next step, summation over /„_i, is performed similarly to (6.6). The only 
difference is in the presence of the additional geometric progression r*"-^. The 
summation yields 

/ |1 - ■ - ' " - 

const — 7 

VI - r \Un-l\ 

Once again, for small variations of Un-i, 1^71-1,71 , w„, this is bounded by a constant 
times a geometric progression with exponent m„_2,n-i and a ratio that is close to 
1. Induction on n and the presence of in 'Kzz'^^{li) complete the proof of the 
uniform convergence of the series. 

Let us prove (ii). The first step is the same: the sum (6.2) is split into the outer 
sum over nonnegative integers h,. . . ,ln and mi2, . . . , m„_i,„ satisfying (6.4), and 
the inner sum over Young diagrams A(i), fi{i, i + I) with |A(i)| — li, \fi{i, i + 1)| = 
mj,i+i restricted by A(i) G Vi and /u(i,« + 1) G Cj^j+i. As was mentioned above, 
if all Vs are equal to Y, the inner sum yields (6.5). Since all the summands are 
nonnegative, with arbitrary T^'s the inner sum yields (6.5) multiplied by a constant 
which depends on Vs and m's, does not depend on ^ and 77, and is between and 
1. It is also worth noting that these constants are polynomials in z and z' because 
the inner sums are always finite. 

As we replace ^ and tj by £^ and st] in the expression obtained by using the 
formulas of Lemma 4.8 in (6.5), we use the estimates 



(1 - eri)e^J 

1 — ST] 



ZZ +771 



(1 - ei)eri 



II 



{I - m)W. 



< const 



1 — er] 



{c + x) 



{I - m)\ 



< lei' • const 



l—m 
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where the bounds are uniform in e varying in the unit disc |e| < 1. These estimates 
together with |7rc,£{(/)| < |£|' const' imply that (6.5) is bounded by 

|£|ii+-+i„-mi2 m„_i,„ . const'^+-+''' . (6.7) 

Let us show that the radius of convergence of the power series in |e| obtained by 
adding expressions (6.7) with all nonnegative (Zi, . . . , /„, mi2, . . . , m„_i,„) satisfying 
(6.4), is positive. Indeed, it is not hard to see that for any numbers satisfying the 
system of inequalities (6.4), one has 

\-ln- mi2 m„_i,„ > max{Zi, . . . , 

We can split the power series into parts according to which of Zi's is the largest one. 

It suffices to check the convergence of each part. Let us take one of such subseries, 
say, assume that I = Ij = max{/i, . . . , Z„} for some j. Then for |^| < 1 

|^|h+-+i„-mi2— -m^„_,^„ < const'^+-+'" < max(l, const"'). 

Finally, the number of sets of 2(n — 1) nonnegative numbers (not including 
Ij = I) and {mi,i+i} bounded by I from above is {I + l)^^""^). Thus, our series is 
majorized by 

oo 

^(; + l)2(n-l) const"' |£|', 

1=0 

and this series converges for small enough |£|. 

Thus, we have verified that the series of expressions (6.5) absolutely and uni- 
formly converges when |e| is small enough. Clearly, the multiplication of the terms 
of this series by constants between and 1 which we mentioned at the beginning 
of the proof of (ii) (these are the inner sums over Young diagrams) does not affect 
the convergence. This proves the first statement of (ii). 

The second statement of (ii) is easy: the terms of (6.5) with ^, tj replaced by 

ST], have a zero of order at least Zi H \-ln — toi2 m„_i^„ at e = 0; it comes 

from the factors (e^; — tr/i-i.!)'''^""*"^^' in P^'s and from in 7rc.^(Zi). Thus, only 
finitely many terms contribute to a fixed Taylor coefficient. Each of these terms 
involve polynomial expressions in z,z' and expressions of the form (1 — e^iY^ , 
(1 — e7}is+i)~^^ , and their Taylor coefficients at e = are also polynomials in 



7. Dynamical correlation functions 

7.1. Definitions. Consider a continuous time stochastic process A{t) with the 
state space Y (all Young diagrams). As in §3 we view the Young diagrams as point 
configurations (=subsets) of Z' = Z + ^ via 

A e Yh^ X(A) = (a;i,a;2,...) C Z', = Aj - i + ^ , i = l,2,... 

Then A{t) is equivalent to the corresponding process with values in point configu- 
rations in Z'; let us denote this process by X_{t). 
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For any n = 1, 2, . . . define the nth dynamical correlation function of n pairwise 
distinct arguments {ti,Xi), . . . (t„,a;„) e {t„ii„,tmax) x Z' by 

Pn{ti, xi;t2,X2; ■ ■ tn, Xn) = Prob {X_ (ti) contains Xi for alH = 1, . . . , n} . 

In other words, the dynamical correlation fimctions describe probabilities of 
events of the following type: for given time moments si < • • • < Sm and given 
finite sets Yi,. . . ,Ym, the random point configurations X_{si), . . . ,2^{sm) contain 
Yi, . . . , Ym, respectively. Thus, the notion of the dynamical correlation functions is 
a hybrid of the finite-dimensional distributions of a stochastic process and standard 
correlation functions of probability measures on point configurations. 

Clearly, the dynamical correlation functions uniquely determine the finite di- 
mensional distributions of the process and, thus, the process itself. The reason 
why we are interested in these quantities is the same as in the "static" (fixed time) 
case: As we take scaling limits of our processes, for the limiting object the notion 
of the weight of a point configuration does not make sense anymore. Thus, the 
probabilities of the form Prob{X (si) = ^i, • • • ,X_{sm) = X_m} do not have any 
meaning in the limit while the scaling limits of the correlation functions are well- 
defined and, moreover, carry a lot of useful information about the limit process. 

We say that the process X_ (t) is determinantal (cf. §3) if the exists a kernel 

^ ' {^{tmin tutax^ X Z ) X (^(tijtin j ^max^ X Z ) ^ C 

such that for any n — 1,2, . . . 

Pn {tl-^Xx, . . . , Xn , ^n) — dct [^7^(tj , Xi , tj , Xj )j ^ j — \' 

As in the "static" case, if such a kernel exists then it is not unique. In particular, 
transformations of the form 

K{s, x; t, y) — > ^-j^ K{s, x; t, y) (7.1) 

do not change the correlation functions. 
7.2. Main results. 

Theorem 7.1. Let {z, z') he a pair of admissible parameters and ^(•) he an ad- 
missible curve. Consider the Markov process A2,^j,',€(-) defined in %4i O'nd denote by 
^2,z',{( ) the corresponding process with values in the space of point configurations 
in 11 . Then the process X_z,z',^{-) is determinantal. 

Recall that in (2.1) we introduced the functions a{x; z , z' , which form, for 
any ^ e (0, 1), an orthonormal basis in These functions were defined under 

the condition that {z, z') belong to either principal or complementary series. 

Theorem 7.2. Assume that {z, z') is either in principal or complementary series. 
Then the kernel 

Kz,.',ii.){s,x;t,y) = ±J2 e^''^'-'H±a{x;z,z',^{s))i,±a{y;z,z',at)) (7-2) 
aez^ 

with "+ " taken for s > t and "— " taken for s < t, is a correlation kernel of the 
process Xz^z',i(-)- 
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Theorem 7.3. The correlation kernel (7.2) can also be written in the form 

Kz,z',i{-){s,x;t,y) = e5(''-'V2,0'(2;,2/)^2,z',c(-)('S.a;;t,2/) (7.3) 
where, as in (2.6), 

^T{x + z + i)r(x + z' + i)r(t/ + z + i)r(y + z' + ^) 



T{x + z' + \)T{y + z+\) 



and the kernel K^z,z',({-){s,x;t,y) can be written as a double contour integral {set 



Kz,z',i{-){.s,x;t,y) 

_ V(i-C)(i-^) 



{2mf 



LJ^ LJ. 



-y-2 



2 



■ dojiduj 



2 



* (wi - C) (w2 - r?) - (1 - C^i) (1 - r?W2) 

(7.4) 

with the contours {wi} and {^2} o/wi and u)2 satisfying the following conditions: 

• {uJi} goes around in positive direction and passes between C, and C,^^ ; 

• {(JJ2} goes around in positive direction and passes between rj and r]~^; 

• if s>t then the image of {wi} under the fractional-linear map 

uj {e^-^Tj - C) + 1 - e^-^Cr] 



(j{e^ * — C??) + 7? — *C 



(7.5) 



is contained inside {102}; 

• if s < t then the domain bounded by {0^2 } does not intersect the image of {oji} 
under the map above. 

The kernels I£z,z',(,{-){s,x;t,y) and I£z,z',(,{-){s,x;t,y) are equivalent. Namely, 
they are related by a "gauge transformation" (7.1), 

Kz,z',i{-){s,x;t,y) = {T'"" \^'^} Kz,z',i{-){x,y), x,y e Z', 

where 

fz,z\s,x)= = (7.6) 

^r{x + z+^)T{x + z' + l) 

The kernel K_^^^^^(^.^{s,x;t,y) can serve as a correlation kernel for all admissible 
values of parameters {z,z'), including the degenerate series. 

Comments. 1. The fractional- linear transformation (7.5) arises from the condition 
that the denominator in the integral representation (7.4) has to be nonzero. Solving 
the equation denominator=0 with respect to ui yields the right-hand side of (7.5) 
with u> = UJ2. 
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2. It is not a priori clear why the needed contours {oJi\ and {^2} exist. Let 
us show that it is indeed so. Set q = e^~*. Note that (7.5) maps ( r]~^ and 

Consider the case q> 1 (i.e. s > t) first. Let us take a circle with center at the 
origin and radius r slightly smaller than as {u)i}. Then its image is again a 
circle which is symmetric with respect to the real axis and which passes through 
the images of r and —r. The image of r is close to the image of which is rj, and 
the image of — r is close to the image of — which is equal to 

-g^(C + C-^) + 2 

-9(C + C-^) + 2r?- 

Since C and 77 are strictly between and 1, we immediately see that the denominator 
is negative, and the whole expression is < 77. Thus, the image of {coi} is a finite 
circle that lies to the left of 77 plus a small number. Clearly, there exists {0J2} that 
passes between rj and rj~^ and encircles both and the image of {wi}. 

Let us consider the case q < 1 now. As {coi} we again take a circle with center 
at the origin but with radius slightly greater than (. Then its image is a circle 
which is symmetric with respect to the real axis and which passes through images 
of points that are close to the image of ( which is r]~^, and to the image of — C 
which is 

c+r'-2gr? 

rj{C + C-^)-2q- 

If the denominator is negative then the whole expression is negative, and there exists 
a contour {^2} inside this circle that passes between rj and r]~^ and goes around 
the origin. If the denominator is positive then the whole expression is > r]~^, and 
{C02} can be a circle of radius between r/ and ri~^ with center at the origin. □ 

Theorems 7.1, 7.2, and 7.3 are generalizations of Theorems 3.1, 3.2, and 3.3, 
respectively. 

7.3. Proof of Theorems 7.1, 7.2, and 7.3. The ideas used in the proof are 

similar to those of §3. The first step is to consider the case z = N G {1,2, . . .}, 
z' = N + a, a > —1. Then the state space of our Markov process ^z,z',(, is smaller 
than the whole set Y. Namely, since MArjv+Q,4(-^) vanishes if (the number of 
nonzero rows of A) is greater than N, our process lives on the set Y{N) of Young 
diagrams with no more than A'^ rows. Consider an embedding of Y(A'') into the set 
of iV-point subsets of Z_|_ given by 

XeY{N)^ X{X) = (xi,...,xn), x, = X, + N-i, i = l,...,N, 

and denote by ^Ar.jv+a,? the corresponding stochastic process on the space of Ap- 
point configurations in Z+. Obviously, the processes ^jv,jv+a,{ and XN,N+a,^ are 
equivalent: if A e Y(A'') then X_{X) = {xi,X2, ■ ■ ■) and X{X) = {xi, . . . ,xn) are 
related by 

_(x^-N+^, i=l,...,N, 
i>N + l. 

This implies that the dynamical correlation functions pi of :Kjv,JV+a,{ and pi of 
XN,N+a,$ are related by 

Pi{n,xi; . . .■,Ti,xi) = pi{ti,xi; . ..■,ti,xi) 
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with Xi = Xi -\- N — ^ Z+ for i = 1, . . . , Z. 

Define tlie extended Meixner kernel K^^^^^^ by, cf. Lemma 3.4, 



x-tv)-{ *' ^' *' 



witli a;, J/ e and 



^Jf,t£T)'^(«' x;t,y)=Y, e™(«-*)9Jl™(J; a, a, mi (7-7) 

m=0 

CSO 

<Sr(^'^5t,y) = - E e'"(^-*)sat™(J;a,^(s))£at™(y;a,e(t)). (7.8) 



m=N 



Lemma 7.4. T/ie process ^jv,jv+a,{(-) determinantal. Its correlation functions 
have the form 



p; (n , xi ; . . . ; r; , X; ) = det -K'^^^'?/':) (rj ,Xi;Tj,Xj) 



1 = 1,2,. 



We postpone the proof of this lemma till §7.4. 

The next step is to connect the extended Meixner kernel /i'Me™er g^^^ ^j^g kernel 
K of Theorem 7.2 (cf. Lemma 3.7). 

Lemma 7.5. We have 

K^tii]{s,x;t,y) = e(^-5)(-*)^;v,Ar+a,|(.)(s,x;t,y) (7.9) 
with x = x + N- Ig Z+, y = y + N-\& Z+. 

Proof. We argue as in the proof of Lemma 3.7. For s>t take (7.7) and change the 
summation index m = N — a — \. Then, see Proposition 2.8, 

^m(s-t) ^ g(JV-i)(.-t) . ^a{t-s)^ ^) ^ AT, AT + a, Cj 

with X = X + iV - i. Furthermore, Va(a;; N,N + a,i) = Q \{ a = N + \,N +\, . . . 

because of the factor r(z — a + i) in (2.1). This yields (7.9). For s < t the argument 
is similar; it uses (7.8), the summation index change m = N — a — \ and 

em(.-t) ^ g(iv-i)(.-t) . e-a(.-t)^ ^) ^ jv, TV + a, 0. 

□ 

Lemma 7.4 and Lemma 7.5 imply that the correlation functions pi of X;v,Ar+c«,^(-) 
can be written as 

Pi(Ti,xi; . . . ;t;,x;) = det ^;v,Jv+a,|(-)(Ti,a;i;Tj,Xj) , Z = l,2, 

(7.10) 

if Xj > -Af + i for alH = 1, . . . , « (indeed, the factor e(^-^)(*-*) in (7.9) does not 
affect the correlation functions). 

The next claim is a counterpart of Lemma 3.9. 
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Lemma 7.6. Assume that 

• either {z, z') is not in the degenerate series and x,y &Z' are arbitrary 

• or z = N = 1,2, ... , z' > N — 1, and both x, y are in Z+ — N + ^. 

Then the kernel K_z^z',${-){s, x;t,y) of Theorem 7.3 is related to the kernel K_z^z',i{-){s, x;t,y) 
of Theorem 7.2 by equality (7.3). Or, that is the same, by the "gauge transforma- 
tion" 

Kz,z',i(-){s,x;t,y) = {f'"" ['^''^l Kz,z',i(-){s,x;t,y), 

jz,z' y) 

where fz.z' is defined in (7.6). 

Proof. Applying formula (2.4) of Proposition 2.3 to and setting a = A; + ^ we 
obtain 



^l:k+dx;z,z',£,{s)) 



ijk+h{y;z,z',i{t)) 



r{x + z + i)r(a; + z' + i) y r(z' - k){i - ^(s)) 



r(z - k)T{z' - k) 

-z' + k 



T{x + z' + \] 



i-^j{s)uj^) ' (i-vi(^^r' 



r{x + z + i)r(x + z' + l)y T{z- fc)(i - e(t)) 



r{z - k)r{z' - k) 

-z+k 



X i (i - ^^^^2) (1 - Vm^2') 



T{x + z + \) 
^'-'^-i -y-fe-i dW2 



U)2 



(As in the proof of Proposition 2.3, we used the symmetry ^a{x; z, z' , ^) = 'il>a{x; z' , z, ^) 
to get the second relation.) 



Note that both sides of (7.3) are real analytic fimctions of q 



on g > 1 



and Q < q < 1 with all other parameters [z, z' ,S,{s),£,{t),x,y) being fixed. Thus, 
it suffices to prove (7.3) for q large enough in the case of g > 1, and for q small 
enough in the case q <1. 

Let us consider the case q> 1. Substituting the integral representations above 
in (7.2), we observe that the summation index fc = a — ^ € Z+ enters the resulting 
expression only as the exponent in 

q {uj, - vM) {^2 - vew) 

For large enough q the geometric progression X^j,>q with this ratio converges uni- 
formly on any fixed contours {ci;i},{a;2}. Computing the sum yields (7.3), (7.4). 

Similarly, in the case g < 1 the computation reduces to summing the geometric 
progression with the ratio 

q (a;i - V^M) (^^2 - ^/m) 
(1 - ^)u^i) (1 - VM'^2) 

which always makes sense for small enough q. 

For large q the image of any finite contour under (7.5) is concentrated near 
rj = -\/^(t), and for small q such image is concentrated near r]~^. These points 
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are always inside/outside of any contour {0^2} provided by Proposition 2.3. The 
conditions on contours in the statement of Theorem 7.3 reflect the deformations 
of the contours of Proposition 2.3. Note that the denominator of (7.4) must stay 
away from zero while deforming the contours, which means that the image of {wi} 
under (7.5) is not allowed to intersect lo2- □ 

The relations (7.10) and (7.3) yield the following formula for the correlation 
functions of XN,N+a,$,'- 

Pi{ti,xi;...;ti,xi) =det K;,^^,^^(^.){Ti,Xi;Tj,Xj) , 1 = 1,2,..., (7.11) 

L -l«,J = l 

where (z, z') = {N, N + a) and Xi + TV — ^ g Z+ for alH = 1, . . . , Z. In order to 
prove Tlicjorcm 7.3, wc need to cxtcncl this formula to arbitrary admissible {z, z') 
and arbitrary Xi G Z'. We will do that by means of Proposition 6.2. 

Up till now the time moments n , . . . , r; were arbitrary, they were not ordered 
and some of them were allowed to coincide. Let us denote by ti, . . . . i„, n < I, the 
same numbers but ordered and without repetitions. Thus, each is equal to one 
and only one tj. As in §6, we set 



?i = ^iU), I <i<n, rji^^+i = e*- ^/Mi+i- 

Then in the notation of §6, see (6.2) and below, pi{t\,Xi; . . .■,Ti,xi) is equal to 
F(^, 77) with a suitable choice of the sets 2?i , . . . , I>„ and Vx2, ■ ■ ■ , 2^n-i,n- Namely, 
we take all equal to Y, and the set Vi is determined according to the following 

recipe: take all numbers j G {1,. . . ,1} such that tj = U, then the corresponding 
points Xj G Z' must be pairwise distinct. Then 

I'i = {a G Y I X (A) contains all Xj such that Tj — ti]. 

Proposition 6.2 says that pi{ti,xi\ . . . :ti,xi) is a real analytic fmiction in ^.r], 
and after the substitution (^, ry) 1-^ (e^, er/) the corresponding function in e can be 
analytically continued in a neighborhood of £ = 0. Moreover, its Taylor coefficients 
at this point are polynomials in z and z' . 

Now let us look at the right hand side of (7.11) with N,N + a replaced by z, z' . 
The values K_z^z'4{-){TijXi;Tj,Xj) of the kernel are given by (7.4). This formula 
involves y^S^Ti), \/£.{Tj), and e^^~'^^ , which is expressible through v^'s and ?7i_i+i's 
in a polynomial fashion. Moreover, e^^~'^^ do not change if we scale ^ and r]hy e. 

The integral representation (7.4) implies that K_ z,z' ,^{-)iTi, xt; tj , xj) are real an- 
alytic fimctions in v^'s and jjj^i+i's. Further, if we scale ^ and rj by e, then (7.4) 
viewed as a function in S = , extends to an analytic function in a small enough 
disc {6 : \d\ < const}. ^ Moreover, its Taylor coefficients at (5 = are polynomials in 
z and z' because the Taylor coefficients of (1 — u)" at u = are polynomials in k. 

We conclude that both sides of (7.11) are uniquely determined by their values 
for {z, z') = {N, N + a) with any large enough N and any a > —1. This completes 
the proof of Theorem 7.3. 

Theorem 7.2 is a direct corollary of Theorem 7.3 and Lemma 7.6. Theorem 7.1 
follows from Theorem 7.3. □ 



'^It is worth noting tiiat for small e we can choose the contours of integration in (7.4) which 
would be independent of e; it suffices to consider suitable circles centered at the origin. 
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7.4. Eynard Mehta theorem and the proof of Lemma 7.4. Here we state 
the Eynard-Mehta theorem [EM] in a form which is convenient for our purposes 
and show that Lemma 7.4 is a corohary of this theorem. 

Let TO be a fixed natural number and let the index k range over {!,..., m}. 
Consider the Hilbert space taken with respect to the counting measure on 

the set Z_|_ = {0, 1, 2, ... }. Assume that for each k we are given an orthonormal 
basis {^!'fe,Ti}n=o,i,... of real-valued functions in ^^(Z_|_). Next, assume that for each 
k = 1, . . . , m — 1 and each n = 0, 1, . . . we Rre ci given Ei number cj^^k-\-i',n 

> 0. 

As n — > c», these numbers have to decay fast enough to make convergent certain 
infinite sums specified below. Finally, we will impose on these data certain positivity 
conditions, see below. 

We aim to construct a probability measure on collections (Xi, . . . ,Xjn), where 
each Xk is an arbitrary A^-point subset in Z_|_ and A' is a fixed natural number. 
This measure can be regarded as a Markov process with "time" k = 1, . . . , m, the 
state space being the set of AT-point subsets in Z+ . 

The construction goes as follows. For an arbitrary A^-point subset X = {x\ < 
■ ■ ■ < Xn) C Z+ we introduce the N x N matrix (/)fe(X) with the entries (j)k,i{xj), 
where the row index i takes values in {0, . . . , A^ — 1}, and the column index j takes 
values in {1, ... , A'}. 

As X ranges over all A'-point subsets of Z_|_ , one has 

^(det.^fc(X))2 = L 

X 

The proof follows from the Cauchy-Binet identity and orthonormality of (j}k n's, cf 
(5.8). 

Thus, for each fc = 1, . . . , m we have a probability measure dk on A'-point subsets 

in Z+ which assigns to a subset X its weight (det (/)fc(A'))^. The measures ak are 
the 1-dimensional distributions for our future Markov process. 
For each A; = 1, . . . , m — 1 we set 

oo 

Vk,k+l{x,y) = '^Ck,k+l;n<Pk,n{x)(f)k+l,n{y), t/ € Z+ , 

where the sum is assumed to be convergent. Since {(j)k,n} is an orthonormal basis 
for each fixed k, we have 



(l)k,n{x)Vk,k+lix,y) = Ck,k+l;n<f>k,niy), 



'+ 



(7.12) 

Vk,k+l{x,y)(l>k+l,n{y) = Ck,k+l;n(pk,n{x). 

For arbitrary subsets X = {xi < ■ ■ ■ < Xn) and Y = (jji < • • • < y^) we form an 
N X N matrix Vk,k+i{X,Y) with entries Vk,k+i{xi,yj), and we set 

,^ ^, det0fe(A)detwfc,fc+i(A,y)det</)fe+i(y) 
(Tk,k+i{X,Y) = 

n Cfe,fe+l;»i 
n=0 
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Once again, using the Cauchy-Binet identity and (7.12), it is not hard to show 
that for any fc = 1, . . . , m — 1 one has 

'^(^k,k+i{x,Y) = (Tk{x), ^(7fc,fe+i(x,y) = (Tfe+i(y). 

Y X 

Assume that ak^k+i{X,Y) > for all X and Y. Then we may regard ak^k+i 
as a probability measure on couples {X, Y) with marginal measures ak and ak+i ■ 
This is the "one step" 2 dimensional distribution of our Markov process. Wc also 
assume that det(j)k{X) does not vanish. Then we define the "one step" transition 
probability function as follows 

t> fl.l.,^^ q-fe,fc+i(X,y) det Vk,k+i {X, Y) det (j)k+i (Y) 

Px,Y{k,k + l)= ^^^^^ = — . (7.13) 

det4>k{X) n Cfe,fe+i;„ 

Wc regard this as a matrix P{k,k + 1) whose rows and columns are labelled by 
A^-point subsets. 

Finally, we define a Markov process X{k), where the "time" k takes values from 
1 to m and X{k) is an point subset of Z+, using the initial distribution cti and 
the "one step" transition probabilities (7.13): 

Prob(A'(l), . . . , X(m)) = (Ti{Xi)Px(i),xi2){^, 2) . . .Pxim-i),x{m){m - 1, m). 
For arbitrary indices k, I such that 1 < A; < / < m we set 

Ck,l;n = Ck,k+l;n Cfc+l,fc+2;n • • • Cl-l,l;n , Tl = 0, 1, . . . , 

oo 

Vk,i{x,y) = 

n=0 

Theorem 7.7 (Eynard— Mehta [EM]). Under the above assumptions, let us re- 
gard the Markov process X{k) as a probability measure on mN -point configurations 
X = {X{1), . . . ,X{m)) in the space {1, . . . ,m} x Z_|_. Then this measure is deter- 
minantal, and its correlation kernel has the form 

N-l ^ 

K{k,x;l,y)= ^ (j)k,i{x)(j}i^i{y), k>l, 

[where we agree that Ck,k;i — 1) '^"'^ 

N-l 

K{k,x;l,y) = Ck,i;i(l)k,i{x)(l)i,i{y) - Vk,i{x,y) 

i=0 
oo 

= - X] ^'''^•i ^k,i{x)(l)i,i{y), k<l. 

In other words, for any n = 1, 2, . . . we have 

Pn{ki,xi; . . . ; kn,Xn) ■ = Prob{X(fci) 9 Xi for each i = 1, . . . ,n} 
= det [K[ki, Xi; kj, Xj)]^ j_-^ , 
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where {ki,Xi) € {1, . . . , m} x Z and {ki,Xj) ^ {kj , xj) for i ^ j. 
Proof. See [EM], [NF], [Jo3], [TW], [BR]. □ 

Proof of Lemma 7.4- The process Atv jv+q,,^ restricted to any finite sequence of time 
moments ti < ■ ■ ■ < tm fits into this formahsm perfectly. Indeed, we take 

(^k,n{^.) = OT„(.t; a, e(tfc)), Ck,k+i;n = e''^'''-'''+'\ 

Then (5.7) or Lemma 3.4 imply that Mjv,Ar+a,j(t^) is exactly Ufc, and the transition 
matrix (5.2) coincides with (7.13).^ Lemma 7.4 is thus a direct corollary of Theorem 
7.7. □ 

7.5. An interpretation via nonintersecting paths. In this section we interpret 
the stationary process Ajv,jv+a,j in terms of N nonintersecting trajectories of in- 
dependent birth-death processes. This is done using formulas of Karlin-McGregor 

[KMG3] and an idea of Johansson [Jo2]. 

Instead of dealing with AArjv+a,{ we will use the associated process XN,N+a,i 
introduced in the beginning of §7.3. Recall that its state space consists of AT-point 
configurations in Z+. 

In the special case A^ = 1 the process Xi i+q .c is just the birth-death process 
Ni^a,^- We aim to construct A"Arjv+Q.^ directly in terms of A^i+a,^. 

Let us take a large T > and consider a new process YAf.a.^.T introduced as 
follows. This process is defined on the time interval [— T, T]. Let us take A^ in- 
dependent birth-death processes which start at the moment — T at the points 
ai < ■ ■ ■ < On and end up at the moment T at the points 6i < • • • < 6]v con- 
ditioned on the event that the trajectories Xi{t) do not intersect on [—T,T]: 

Xi{t) <X2it) <■■■ <XN{t), -T<t<T. 

The boundary conditions {ui} and {bi} are arbitrary but fixed while the parameter 
T will vary. 

Theorem 7.8. In the above notation, as T ^ oo the processes ljv,a,{,T converge 
to XN,N+a,i sense of convergence of the finite dimensional distributions. 

Proof. Let us fix arbitrary time moments ti < ■ ■ ■ < tk inside {—T,T). Then 
by [KMG3] and Theorem 5.1, the corresponding finite-dimensional distribution of 

YN,a,i,T has the form 

Proh{YN,o.,i,T{t^) = iVi^ <---<yN) for alH = 1, . . . , fc} = 

det[v.T,tM,,y^'^)] det[vt,,M'\yf)] ' ■ ■ det[vu_,,tdyt'\yf)] det[^t.,r(yf \ 6,)] 

det[v-T,T{ai,bj)] 

where Vs,t{x,y) is given by (5.1) with ^( - ) = ^. 

On the other hand, the finite-dimensional distribution of XN,N+a,i given by 

Pvoh{XN.N+a.dti) = (yP < • • • < 2/^^) for alH = 1, . . . , fc} 
= e^^--^ ^^^^1^ det[.,.a.«, ^)] • • • det[..,_,,,(,f-), ,f ))] 



^Recall that the Young diagrams with no more than N rows are identified with iV-point subsets 
of Z+ via A (Ai -t- iV - 1, A2 + W - 2, . . . , Aat). 
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with (j)n{x) = 9Jt„(x; a, sec Theorem 5.1. 

Note that we have the foUowing asymptotic relation: for arbitrary x'i,x'{, . . . . 



det 



Y^e^4>n{x[)4>n{x';] 



.n=0 



= det[<^^_i^)] det[cPi.,{x'!)] + O (e^^+i) 



as e ^ 0, of. (5.4). Applying this asymptotic relation to V-T,ti, '^tk,T, V-t,t, we 
obtain 

det[u_T,ti(ai,2/} ')] ~ e 5 det[0i_i(aj)] det[0i_i(j/} ')], 

detK„T(2/f \ ~ e- det[</.,_i(yf )] det[<^i_i(6,)], 
det[?;_T,T(ai,?>i)] ~ e-^^'^-^) det[.^i_i(a,)] det[.^i_i(6,)], 

as T ^ +00. This completes the proof. □ 

8. Particle-hole involution 

For any set X and its subset 2) one can define an involution on point configura- 
tions XcXbyXi-^XA2)- This map leaves intact the "particles" of X outside of 
2), and inside 2) it picks the "holes" (points of 2) free of particles). This involution 
is called the particle-hole involution on 2). 

The goal of this section is to give a different description of the z-measures using 
a new identification of Young diagrams and point configurations on Z'. Instead of 
using the configurations 

X{X) = {Xi-i + l\i = l,2,...} 
we will use the configurations 

X{X) = {X (A) n Z'+) U (Z'_ \ X (A)) =X{X)A Z'_ 

which are obtained from X (A) by applying the particle-hole involution on Z'_. 

The parametrization of Young diagrams A by configurations X{X) corresponds 
to considering the Frobenius coordinates of A, see [BOO, §1.2] for details. The 
reason for passing to X{X) is very simple: in the continuous limit y 1 which will 
be considered below in §9, the point process generated by ^(A) does not survive, 
while the process corresponding to X{X) has a well defined limit. 

Observe that X(A') = -X{X) for any A e Y. 

Given an arbitrary kernel K{x, y) on 3C x X, and a subset 2) of 3C, we assign to 
it another kernel, 

K{x,y), a; ^2), 

5xy-K{x,y), a;e2), 

where S^y is the Kronecker symbol. Slightly more generally, given an arbitrary map 
£ : X — *■ M* , we set 

K°'%x,y)=e{x)K°{x,y)e{y)-\ 



K°{x,y) = 
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Proposition 8.1. Let P he a probability measure in point configurations on a 
discrete space X and let P° be the image of P under the particle-hole involution on 
2) C X. Assume that the correlation functions of P have determinantal form with 
a certain kernel K{x,y), 

Pm{xi,...,Xm\ P)= det [K{xi,Xj)], m = l,2,.... 

Then the correlation functions of the measure P° also have a similar determinantal 
form, with the kernel K°{x,y) as defined above or, equally well, with the kernel 
K°'^{x,y), where the map £ : X — > M* may be chosen arbitrarily, 

pm{xi,...,Xm\P°)= dot [K°{Xi,Xj)]= dct [K°'%Xi, Xj)], 

rn= 1,2, ... . 

Proof. The factor £( • ) does not affect the values of determinants in right-hand side 
of the above formula, so that we may take s{-) = 1. Then the result is obtained 
by applying the inclusion/exclusion principle, see Proposition A. 8 in [BOO]. □ 

Later on wc choose the function e( • ) in a specific way (see (8.3) below) which is 
appropriate for the limit transition of §9. 

The main result of this section is a determinantal formula for the dynamical cor- 
relation functions of Az,z'4 computed in terms of X{X). For any n = 1,2,... 
define the nth dynamical correlation function of n pairwise distinct arguments 

(f 1 , Xi ) , . . . (tn , Xn ) € itmin i tmax ) X Y by 
Pn{tl,Xi;t2,X2; ■■ ■■,tn, Xn) 

= Prob {-'^(A) at the moment ti contains x, for alH = 1, . . . , n} . 

Here and in what follows we denote by J^z,z',s, the stationary Markov process 
corresponding to the constant curve ^{t) = ^, where ^ G (0, 1) is fixed. 

Theorem 8.2. Let {z, z') be a pair of admissible parameters not from the degen- 
erate series. Consider the Markov process h^^.z'.^, and denote by Xz^z' the process 
with values in the space of point configurations in 11 , obtained from l^z,z' ,i via 
X^X{X). 

Then the process Xz.z'.^ is determinantal and its correlation kernel has the form 
(in all the formulas below x and y are positive, in ± and =F the upper sign corre- 
sponds to the case s>t, and the lower sign corresponds to s <t) 

Kz,z',i{s,x;t,y) =± ^ e-''^''-*^i}±a{x;z,z',^)i}±a{y;z,z',^), 
aez^ 

Kz,z',i{s, x; t, -y) = ±J2 (-l)^''-*e-»l«-*l V±a(x; z, z', OMv, 0, 
Kz,z',i{s, -x; t, y) = T E (-l)^""'e-»l«-*l V^aCa^; -z, -z', ^)^±a{y; z, z', 0, 
Kz,z',ds,-x;t,-y) = T Yl e'^^'-'^^^aix; -z, -z' ,^)^p^aiy; -z, -z' ,^), 

(8.1) 
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where the fourth formula is valid for s^t, and for s = t we have 

-x; s, -y) = ^ i>a{x; -z, -z', ^)tpa{y; -z, -z', ^). (8.2) 

Com,m,ents. 1. For s = t this kernel coincides with the hypcrgcomctric kernel 
derived in [B02], see also [B04], [B05]. In those papers the kernel was written in 
another, so-called "integrable form" , see Remark 8.3 below. 
2. The kernel has the following symmetries {x,y G Z'): 

K,,,,,^{s,x;t,y) = {-ir^'^^-^^^yK,,,,,^is,y;t,x), 

K-z.^z',^{t,-x;s,-y), Sy^t, 

-x; s, -y), s = t. 



Kz,z',ds,x;t,y) = | 



Proof. We use Proposition 8.1. As the initial kernel we take the expression for 
K_z.z',c(-) given in Theorem 7.2, the set X is the union of finitely many copies of Z' 
which correspond to times at which we evaluate the dynamical correlation function, 
and 2) is the union of the same number of copies of Z'_ . On each copy of Z' the 
function s{-) is chosen in the following way: 

. X r 1' x>0, , ^ 

e(^)= / ^ J (8-3) 

L (-1) a; < 0. 

The statement then follows from Proposition 2.7. The last formula (for s = t) arises 
from the relation 

^=^,y- m '^-aix;-z,-z',^)'ip-a{y;-z,-z',^) 
aez^ 

= 'lpa{x\-Z,-z',£,)lj)a{y;-Z,-z',C), 
aeZ'_^ 

which follows from the fact that form an orthonormal basis, see Proposition 
2.4. □ 

Remark 8.3. Denote by K^.z' ,^{x,y) the kernel that is obtained from the kernel 
I£z,z',i{x, y) of §3 by the procedure described above. That is, Kz^z',^ = (i^z,z',{)°'^ 
with £ given by (8.3). This is a correlation kernel for the z-measure Mz.z'.(^^ cor- 
responding to the map A ^ X{X). Clearly, Kz^z' ,s,{x,y) coincides with the special- 
ization of the kernel of Theorem 8.2 at s = t. Let us abbreviate 

^±i(a;) =i)j^i{x;z,z',^), i^±i{x) = tp±i{x;-z,-z',(,). 

We have for x,y G Z^ (cf. (3.13)) 

^iTI V'- i {x)ip^ {y) - V'l {x)il^- 1 (y) 



z',dx,y) = 



K 



■■,z',dx,-y) = Yzr^ 



x-y 

V- i {x)^- i (y) + V-i {x)^^ {y) 



x + y 



ipi{x)ipiiy)+ilj_i{x)ip_i(y) 

Kz,z',d x,y)- 

, y^zT^ '^-iix)'ip^iy) -■^Piix)^-iiy) 

Kz,z',d-x,-y) = — ^ 2 2 2 

x-y 



(8.4) 
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Indeed, the first three formulas are easily obtained from (3.13), and for the last 
formula we use the symmetry relation 

Kz,z',e.{-x^-y) = K-z-z',d^,y), x,y G z+ . 
These four formulas coincide with the expressions obtained in [B02, Theorem 3.3]. 

9. Limit transition to the Whittaker kernel 

In this section we compute the scaling limit of the kernel Kz,z',i of §8 as ^ y 1 
and the arguments x and y are scaled by (1 — ^). In this limit the lattice Z' turns 
into the punctured real line R* = K \ {0}. 

Let us introduce the continuous analogs of the functions ipa- These new functions 
Waiu; z, z') are indexed by a € Z' and the argument u varies in R>o. They are 
expressed through the classical Whittaker functions Wk,ix{u), see [Er, ch. 6] for the 
definition, as follows: 

Wa{u;z,z') = {V{z-a+\)T{z' -a+\)y'^ u'^ W^_^^^{u). (9.1) 

Since Wre,^(u) = Wk„-h{u), this expression is symmetric with respect to z z'. 
It will be convenient for us to use the following integral representation of Wa- 

r(^' - a + i)e'''(^'-'')u^ 
2n {T{z-a+^)T{z'-a + ^))'^ 

X / C"'''"^""^(l + C)^"""^e-"(f+5)dC. (9.2) 

J +00 

The (standard) notation for the contour of integration means that we start at +oo, 
go along the real axis, then around the origin in the clockwise direction, and back 
to +0O along the real axis. On the last part of the contour we choose the principal 
branch of (^"^+""5, which uniquely determines the values of this function on the 
whole contour. 

This formula is easily seen to be equivalent to one of the classical integral repre- 
sentations for the confluent hypergeometric function ^, see [Er, 6.11.2(9)]. 

Proposition 9.1. If ^ y' 1 and x G goes to +oo so that (1 — ^)x — > u > 0, 

then 

ipa{x;z,z',^) ~ (1 - ^)^Wa{u;z,z'). 



Wa(u;z,z) 



Proof. This statement can be proved in a number of ways, see e.g. [Er, 6.8(1)]. 
We will give an argument which uses the integral representations of ipa and Wa ■ A 
similar argument will also be employed in the proof of Theorem 9.2 below. 

We start with the integral representation (2.4) for tpa- Let us choose as {iv} 
the following contour C(i?, r, ^), where r > is small enough (smaller than the 
distance between ^ and 1 / ^) and i? > is big enough (greater than 1 / ^ + r) : 
The contour starts at the point co = R, goes along the full circle \ui\ = R in the 
positive direction, then along the real line until the point u) = 1/ + r, further 
along the full circle jo; — = r in the negative direction, and back along the 
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real line to u = R. Thus, C{R, r, ^) consists of a "big circle" of radius R, a "small 
circle" of radius r, and a "bridge" between them. 

We now fix R, pick r of order (1 — and take the limit ^ /" 1 of the integral. 
The integration over the "big circle" |w| = i? converges to zero exponentially in 
(1 — ^)~^ thanks to the factor w"^"'*. To take care of the rest of the integral, we 
make the change of the integration variable 

a; = l/VI+(l-OC- 

Then we have 

Note that the second factors in the formulas above are asymptotically equal to 1 
for ^ close to 1 and ( bounded, and are uniformly bounded away from and oo for 
^ close to 1 and ( corresponding to arbitrary u) on the contour. Hence, the rest of 
the integral is asymptotically equal to the absolutely convergent integral 



2TTi 



f (_<^)-'+-^(l+, 
J+oo 



Taking into account the convention stated in Comment 1 to Lemma 2.2 one can 
check that arg(— ^) = —wi on the last part of the contour. Therefore, changing 
(— +'^-2 to +"-2 produces the factor 

^■Ki{z'-a+\) _ ^^ni{z'-a) _ 

Finally, the prefactor in (2.4) asymptotically equals 

(1-0 

{T{z-a + ^)T{z'-a + ^)y 
Thus, (2.4) asymptotically equals 

r(z'-a+i)e'''(^'-'^)u^ 



r(/ - g + ^)u- 



(1-0^ 



27r (r(z-a + i)r(z'-a + i))' 

X / C"'''"^""^(l + C)^"''"^e-"('^+5)dC- 

J+oo 

□ 

Theorem 9.2. Consider the extended hyper geometric kernel Kz^z'^^{s,x;t,y) as 
described in Theorem 8.2. Let £, y 1 and assume that x,y ^ oo inside 11 so that 
(1 — S^)x u, (1 — £,)y V, where u,v gR* . 

Then there exists a limit kernel KYzi{s,u;t,v) on R* x M*; 

\im{l - 0~^K,,,,,^{s,x;t,y) = K^z,{s,u;t,v). (9.3) 
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For s t the formulas for the limit kernel are obtained from formulas (8.1) for 
the kernel Kz,z',i by replacing tpa 's with Wa 's and setting ^ = I: for u,v > 

K^^,{s,u;t,v) = ± ^ e-''\'-'Wa{u;z,z')w±a{v;z,z'), 

KZ'{s,u;t,-v) = ±J2 i-l)^''"^e-''\^-'Waiu;z,z')w^4v;-z,-z'), 

aez^ 

K^'is, -u; t,v) = TYl (-l)^'^"^e"'^'^"*'«'Ta(t*; -z, -z')w±a{v; z, z'), 

a6Z^ 

Kf;,,{s,-u;t,-v)=T J2 e-''\'-'^w^a{u;-z,-z')w:^a{v;-z,-z'). 

Comments. 1. The prefactor (1 — ^)~^ in (9.3) is due to rescaling of the state space 
Z'by(l-e). 

2. The reason of the restriction s 7^ t in above formulas is the divergence of the 
series for ifj*^, (s, ?i; ,s, — w) and K^z'is, —u; s,v). The scries for K^z' {s , u; s , v) and 
KY^,{s,—u;s,—v) do converge and give the correct answer. For s = t there exist 
analogs of formulas (8.4): 

, w 1 (u)w 1 (v) — w 1 (u)w i(v) 



U — V 

w_ 1 {u)w_ i{v) +wi (u)u; 1 (v) 



U + V 

wi{u)wi{v) + w_i{u)w_i{v) 

U + V 

w_ 1 (u)«; 1 (w) — wi {u)w_ 1 {v) 
u — V 



(9.4) 



kY,Au,-v) = ^ 

Here we abbreviate 

■"'zbi (''^) = '*'^±i z, z'\ Wj^\(u) = Wj. 1 (m; — z, — -z'). 

Formulas (9.4) can be derived from (8.4) using Proposition 9.1. They were previ- 
ously obtained in [Bl], [B02, §5]. 

3. In accordance with the terminology of these papers (where the kernel (9.4) 
was called the Whittaker kernel) we call the limit kernel K^^iis, u\ t, v) the extended 

Whittaker kernel. 

Proof. We use Proposition 9.1. In order to prove the theorem, we need to justify 
the interchange of the summation and the limit transition in (8.1). To do this it 

suffices to show that the scries converge uniformly in ^. 
We will prove that each of the four expressions 



\ip±a{x; z, z', 0V'±a(2/; z, z' , 0|, |V'=Fa(a;; -z, -z', Ci'^^aiv, -z, -z' , ^)|, (9.5) 
\tl}±a{x\ z, z', $)V'=Fo(2/; -z, -z', 01, |V'=Fa(a;; -z, -z', 0'4'±a{y; z, z' , ^)\, (9.6) 
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is estimated from above by const(?i, — ^) • g'"', where q > 1 can be chosen 
arbitrarily close to 1, and const (Ujt;) does not depend on a and ^. Together with 
the factors e~"l*~*l in (8.1) this ensures the needed uniform convergence. 

Both expressions (9.5) are estimated in the same way, let us handle the first one. 
We apply formula (2.5) of Proposition 2.3 and we get a double contour integral, in 
which we single out the terms involving a; we observe that all together they can be 
written in the form (^"(1^1, ci.-'2; ^))'^, where k := a— 5 ranges over Z_|_, and cJi and u!2 
are the variables of integration. Let us write down precisely the whole expression 
separately for the upper and the lower choice of sign in the subscript ±a: 

V'a(a;;z,z',C)'0a(?y;-z,z',O 
^ (r(x + z + ^)r{x + z' + i)r(y + z + ^)r{y + z' + _ 

T{x + z' + ^)T{y + z + ^,) ^ 

with 

F++(a;i,a;2;0 = F+(a;i; 0i^+(w2, 0, 

where 

F^M='-^- (9.8) 
For ip-a{x; z, z', ^)tp-a{y; z, z', ^) we obtain a similar expression: 

V'-a(a;; z, z', £,)ip^a{y; z, z' , £) 

(1-0 



(r(x + z + \)Y{x + z' + \)T{y + z + \)T{y + z' + i))i 



T{x + z' + \)T{y + z + \) 



-i-^ / / (^^__(c.„c.2;0)'(l-V^'^i)"'"7l- — 

X (1 - vi-X'" fi - - V' -r^^^-r^^ — — 



with 
where 



(9.9) 



F_{u:,^) = (i^+(^,0)-' = 7^^- (9-10) 

1 - V^w 



Now we need a lemma. 

Lemma 9.3. Let F±{lo,^) he defined by (9.8) and (9.10). For any q > 1 there 
exists a contour C±{^,q) which is of the same kind as in the proof of Proposition 
9.1, and such that 

\F±{uj,(,)\<q Va;eC±(^,g). (9.11) 
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Proof of Lemma 9.3. Recall that in the proof of Proposition 9.1 we used a specific 
family {C(i?, r, ^)} of contours. We will show that it is possible to take C±(^, g) = 
C{R,r,^) with an appropriate choice of parameters R and r (the radii of the "big 
circle" and the "small circle" in C(i?, r, ^)). 

Consider first the case of F+. Fix g > 1 and let u be related to ui by the 
equivalent relations 

~ TP i C\/ 1 ~ 1 + 9V^^ 

uj = F+{uj,^)/q = 



quj - qy/l ' quj + y/l' 

To fulfill inequality (9.11) the contour C+(^, g) must be contained in the image of 
the imit disk \uj\ < 1 under the conformal map uj i-^ lo. Let 5'+(^,g) denote the 
image of the circle = 1; S+{£,,q) is the circle that is symmetric with respect to 
the real axis and passes through the real points 

(these are the images of —1 and 1, respectively). Since we are interested in the 
limit transition as ^ 1 we may assume that ^ is so close to 1 that q > 
Then we have 

Observe that the image of the disk ja;| < 1 is the exterior of S+{^, q) (for instance, 
this follows from the fact that the image of is the point 1/vT which is outside 
S+{^,q)). Now we take C+{^,q) = C{R,r,^), where R and r are chosen so that a 
both the "big circle" and the "small circle" in C{R, r, ^) arc in the exterior of the 
circle S+{^,q): the "big circle" surrounds S+{^,q) while the "small circle" lies to 
the right of S+{^,q). This is possible due to inequalities (9.12) and the fact that 
the distance between the points '^^^ ^-^d is of order 1 — 

The case of F_ is handled analogously. We define lj by the equivalent relations 

q - gV^w gV^w + 1 

Instead of the circle S+{^,q) we have another circle, denoted by S-{^,q), which is 
symmetric with respect to the real axis and passes through the points 

q + V^ q-Vl 



We note that 



The contour C_(^, g) must lie in the exterior of S'_(^,g). We again can take 
C_(^,g) = C{R,r,^) with appropriate R and r. But, in contrast to the case of 
F+, now the "small circle" in C{R,r,^) must surround S-{^,q) (because is 
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inside S-{^,q), see the inequahties above). This requirement can be satisfied be- 
cause the diameter of S-{^, q) (the distance between the points g^y^i and ^^^) 
is of order (1 — 

This completes the proof of Lemma 9.3. □ 

We return to the proof of Theorem 9.2. Let us estimate (9.7). The product of 
the prefactors is asymptotically 

(1-^)U5(^-^')W5(^'-^). 

To estimate the integral we take as contours {wi} and {1^2} the contour C+(^, q^) 
as described in Lemma 9.3. According to Lemma 9.3, on the product of these 
contours, |F++(a;i, a;2; 01 ^ 9> whence we get 



^ f F++(a;i,W2;0'=(l-V^'^i) ' (1 



(1 - ^UJ2 



U}2 



^ 2 , , y 2 



{a;ieC+(|,93)} {a;26C+(|,g3)} 



(1 - V^W2 



UJ2 



z-1 



Wi W2 



Arguing as in the proof of Proposition 9.1 we check that the last integral is uniformly 

bounded as ^ 1. This yields for (9.7) the required estimate of the form const (1 — 
^)q^ with arbitrary q > 1. The estimate for (9.9) is obtained in exactly the same 
way, by using the contours {oJi} = {^2} = C-{£,,q'^). 

The quantities in (9.6) are estimated similarly. We leave the details to the reader, 
and only point out two minor differences. First, while writing the double integral 
representation, we do not need to switch parameters z and z' (as in the proof of 
Proposition 2.3) to get rid of the gamma factors containing a. Second, we have to 
use as {wi} and {cl)2} two distinct contours: either (C+(^,g2) and C-{£^,q^)) or 
(C-(e,g^)andC+(C,g^)). □ 

Our next goal is to give an integral representation for the extended Whittaker 
kernel computed in Theorem 9.2. Let us introduce contours C±{q) for g > 1. They 
can be viewed as limits of the images of the contours C±(^, g) as ^ 1 in the new 
variable C, where 

'^ = i/V!+(i-OC- 

The contour Cj^{q) starts at +00, goes along the real axis, circles around in 
the negative direction, and returns to +00 along the real axis. It has to leave on 
its left the point —1 together with the circle which is symmetric with respect to 
the real axis and passes through the points —q/{q — 1) and —q/{q + 1) (this circle 
contains —1). 
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The contour C_(g) also starts at +00, goes along the real axis, circles around 
in the negative direction, and returns to +00 along the real axis. It has to leave on 
its left the point —1, and it has to leave on its right the circle which is symmetric 
with respect to the real axis and passes through the points — + 1) and l/{q—l) 
(this circle contains 0). 

Note that if C e C+(5) then |C/(1 + C)I < 9, and if C e C-{q) then |(1 + C)/CI < I- 

Theorem 9.4. The extended Whittaker kernel K^^,{s,u;t,v) of Theorem 9.2 for 
s has the following integral representation {u > 0,v > 0) : 

KY^,{s,u;t,v) = e'^*(^+^')(u/w)^e5(«-*) 



0- 



j> ^cr^'(i+circ2-'(i+C2) 



2 



(27ri)2 7 7 ^1 ^ ■ .2 V e-*(l+Ci)(l + C2)-CiC2 

+00 +00 

where for s > t both contours {Ci} and {C2} are of the form C+(e^(*~*)), and for 
s < t both contours are of the form, C_ (e 2 

„w I J. \ (sin(7rz) sin(7rz'))3 w 
K"^ ,{s,u;t,-v) = ^ — ' {u v) 2 e2(« 

sm(7rz') 



f ^cr'(i+circ2"''(i+C2r 



(27r)2 7 7^^ ^ ' e-Hl + Ci)C2-Ci(l + C2) 

+00 +00 

where for s > t, the contour {(1} is of the form C+(e5(''~*)) and the contour {(2} is 
of the form C_(e2(''^'^), while for s < t the contour {(^1} is of the form C_(e2(*~*)) 
and the contour {C2} is of the form C+(e2('^'')); 

K^^,{s,-u;t,v) = -Kf^^,{s,v;t, -u); 
K^,, {s, ~u; t, -v) = KZ,-,, {t, u;s,v). 

The contours may be chosen differently by deforming the contours above so that the 

denominators of the integrands do not vanish. 

Proof. We take the integral representation (9.2) for the functions Wa and plug it in 
into the series of Theorem 9.2. Computing the sum of geometric progression under 
the integral yields the formulas above. The contours are chosen in such a way that 
the absolute values of the ratios of geometric progressions involved are less than 
one. 

As in the proof of Proposition 2.3, in the derivation of the first formula we 
switch z and z' in the integral representation of the second factor, which cancels 
the gamma factors involving the summation index. In the derivation of the second 
formula we do not need to do that, the gamma factors disappear thanks to the 
relation r(x)r(l — x) = tt/ sin(7r.T) . □ 

10. Limit transition to the gamma kernel 

In this section we compute the limit of the extended hypergeometric kernel 
Kz,z' ,(,{s,x]t,y) as ^ 1 and scaling of time s = (1 — ^)cr, t = {1 — with 
finite cr, r € M. 
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Theorem 10.1. There exists a limit of the extended hypergeometric kernel 
Kl—{a, x; T, y) = lim ^,,,,,^((1 - x; (1 - ^r, y) 

where x,y G Z' , a, r G R. 

For a > T, the correlation kernel can be written in two different ways: as a 
double contour integral 

Kl~{a,x;r,y) 

r(-z' -x + i)r(-z -y+ i)e-'^i(^+^')(_i):';+a+i 

(r(-z -x + i)r(-^' -x + i)r(-z - y + i)r(-^' -y + h))'' 

0- 0- 



""Wf f f l + (a-r)+Ci+C2 ^ ^ 

+00 +00 

and as a single integral 

r+oo 

Kf^^'^ia, x; T,y)= / e-^^^-^^w,{u; -z, -z')wy{u; -z, -z')du. (10.2) 
Jo 

The values of the kernel for a < t are obtained from the above formulas using the 
symmetry property 

K~{a,x;T,y) = {-lf+y+'K^—,{T,-x;a,-y), a + r. 
Comment. For a = t the kernel _^s'^™™^ has a simpler "integrable" expression 
^ — (x;y) 



-gamma. _ sin(7rz) sin(7r2;') 



TT sin(7r(z — z')) 
{r(z + X + \)V{z' + a; + i)r(z + y + i)r(z' + y + i)}"'^' 

, (10.3) 



r(z + X + i)r(z' + 2/ + i) - r(z' + X + i)r(z + y + i) 



see [B05, Theorem 2.3]. We called this kernel the gamma kernel. The more general 
kernel i^f'^"™''(c7, x; r, y) of Theorem 10.1 will be called the extended gamma kernel. 

Proof. We start with the series representation of the extended hypergeometric ker- 
nel, Theorem 7.2. Using the symmetry relations of Propositions 2.5 and 2.7, we 
rewrite the kernel in the following form 

K;,^^,^(^{s,x;t,y) 

{ e-«l^-*IVx(a;-z,-z',OV'y(a;-^,-^',0, s>t, 

^ e-'^l«-*l^_,(a;z,z',e)^-j/(a;^,2',0, s<t. 
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This formula implies that if we prove the statement for cr > r, then the case a < t 
will follow by symmetry. Thus, we continue with the assumption s>t and therefore 
cr > T. 

Formula (10.2) is the limit variant of the formula for K_z,z' ,(,{s,x;t,y) above. 
Indeed, by virtue of Proposition 9.1, 

'4}x{a.;-z,-z' ,£,)'4}y{a;-z,-z' ,£,) {I - £,)wa:{u;-z,-z')wy{u;-z,-z') 

provided that a ^ (1 — ^)~^u. The factor (1 — ^) is responsible for turning the sum 
into an integral over u; this sum is just an approximation to the integral. 

This empirical argument needs a rigorous justification. It is simpler to turn the 
series representation for K_z,z',i{s,x;t,y) into a double contour integral and then 
pass to the limit in the integral. The limit integral will be identified with the 
right-hand side of (10.2). 

Using formula (2.5) with appropriately changed parameters, we obtain 

T{-z'-x+^)T{-z-y+^) 



{T{-z -x+ ^)T{-z' -x + ^)T{-z -y+ i)r(-^' - y + ^))' 



(l - v/fwa) 



^ _ , 1 

i — I UJ-, U/o 

t02 J Wi W2 



with contours {wi} and {^2} chosen as in Proposition 2.3. Now we want to sum the 
geometric progression inside the integrals. The ratio of the geometric progression 
is e^~^u>^^u)2^ . In order to justify the interchange of summation and integration 
we need to ensure that the absolute value of this ratio, as a function in wi, (^2, is 
bounded from above by a constant strictly less than one. This is easy to arrange 
by requiring, for example, that both contours contain the unit circle. 
Performing the summation, we obtain 

Kz,z',(.{^-,x;t,y) 

eh{^-t)T{-z' -x + \)T{-z -y + \) 



{T{-z -x+ \)T{-z' -x+ \)T{-z - y + \)T{-z' - y + \)) ' 

\ J \ J e^~*UJiLJ2 — 1 

Recall the notation C {R, r, ^) for certain type of contours introduced in the proof 
of Proposition 9.1. We assume that both integration variables range over such a 
contour with R being a fixed number greater than 1, and r being of order 1 — ^, 
and such that l/\/f — r > 1. 
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Let us split each of the contours into two parts: the first one is the "big" circle 
= R, and the second part is the rest. If both ui and u!2 range over their big 
circles then the integrand is uniformly bounded, and the prefactor (1 — sends the 
whole expression to zero as ^ — > 1. 

If one of the variables, say, u!2 ranges over its big circle and coi ranges over 
the second part of its contour, we observe that all the factors of the integrand 
involving uj2 are uniformly bounded. The absolute value of the remaining part of 
the integrand 



, 1 



(i - V^Wi) 
is uniformly bounded by 



const • 



with LOi ~ + (1 — £,)Ci, where wc used the same argument as in the proof 

of Proposition 9.1. Thus, our double integral is bounded in absolute value by the 
following one-dimensional integral in : 



const •(1-0^^''"'^^ 



-\-X- 



1 



Ci 2^i + Ci)-^-"-2dCi 



with _ 

Hence, our expression is bounded by 



const •(! - ^f(^'--) j\f^''~'^~'dCi 



which is either bounded by a constant (if ^{z — z') ^ 0) or by | ln(l — ^)| (if 
3f?(z — z') = 0). In both cases, the prefactor (1 — ^) in the integral representation 
for K_z,z',$,{s, x; t, y) sends the whole expression to zero. 

The only asymptotically significant part of the integral comes from the case when 
both uji and UJ2 vary over the second parts of their contours. Making the change of 
variables 

'^i = l/V^+(l-OCi, a;2 = l/v^+(l-0C2, 

and arguing as in the proof of Proposition 9.1, we conclude, using the asymptotic 
relation 

1-^ 1 

e^-*uJiU2 - 1 " 1 + (cr - r) + Ci + C2 ' 

that the limit value of the kernel is given by the right-hand side of (10.1). Note 
that the integral in (10.1) is absolutely convergent. To see this we use the estimate 

|1 + (ct - r) + Ci + C2I > const -1(111(2 T"'' 

which holds for any Ci)C2 on our contours and any v G (0,1). We apply this 
inequality with z/ = i + lR.[z' — z)/2. The fact that v G (0, 1) follows from our basic 
assumptions on z, z', sec §1. 

Thus, we have proved the integral representation (10.1). To see the equivalence 
of (10.1) and (10.2) we substitute the integral representation (9.2) into (10.2) and 
integrate explicitly over u. □ 
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